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X-RAY ABSORPTION IN CHEMICAL ANALYSIS AND CONTROL 
By Herman A. Liebhafsky 


Research Laboratory, General Electric Company, Schenectady, New York 


Three X-ray processes are potentially useful in analytical chemistry: the 
diffraction of X rays by crystals, the emission of characteristic spectra, and 
the absorption of X rays by matter. The first process can establish struc- 
ture and composition. The others establish composition only. In the 
past, the analytical chemist has used each process to a very limited extent, 
because there was no convenient way to measure the intensity of an X-ray 
beam with sufficient speed and precision. Of the three, he has logically 
found that X-ray diffraction is the most useful for the following reasons: 
first, he had many ways of determining composition, but few of establish- 
ing structure; second, it is often possible by X-ray methods to establish 
structure, but usually it is impossible to do quantitative analysis without 
making a precise intensity measurement. As a consequence, the unique 
contribution that X-ray absorption methods can make to chemical analysis 
and control, though long known to physicists, has not been generally appre- 
ciated. 

Owing to marked and relatively recent progress in measuring the intensity 
of X-ray beams, this situation is now changing. There are available today 
three X-ray detectors that convert X-ray beams into electric currents pro- 
portional in magnitude under the simplest conditions to the intensities of 
the beams. ‘These devices are the ionization chamber, the Geiger counter, 
and the photoelectric X-ray detector. The ionization chamber, the oldest 
of the three, has found little or no recent application in chemical analysis 
or control, nor, in view of the generally greater suitability of the other two 
devices, is it likely to do so. The analytical chemist can proceed on the 
basis that the Geiger counter or the photoelectric detector, both of which 
will be discussed subsequently, are capable of meeting any reasonable re- 
quirements that he is likely to impose. He is thus in a position to realize 
the advantages of X-ray absorptometry, some of which accrue from the 
facts that it can be accomplished without noticeably altering the sample, 
and that a single measurement can often be made in a matter of seconds 
once the sample is in the beam. 

The advantages and limitations of X-ray absorptometry are in general 
deducible from the known characteristics of X rays. The reader is referred 
to several excellent books'~ for a thorough discussion of these. This brief 
review, being largely qualitative, will present only the necessary minimum 
of fundamental information about X rays, and it cannot discuss all the 
advantages or limitations of X-ray absorptometry. (For a more complete 
treatment, see bibliographical reference 5 and others listed in the bibliog- 
raphy therein.) 

X-Ray Absorption Described. The outstanding characteristic of X rays, 
from which certain others derive, is their high energy, or small wavelength 
(near 1 A). Because of this characteristic, the absorption of X rays usually 
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involves only electrons near the atomic nucleus. The absorption of ultra- 
violet, visible, or infrared rays, on the other hand, affects the electrons that 
determine the chemical properties of the elements. X-ray absorption is 
thus essentially an atomic process and virtually independent, therefore, of 
the chemical or physical state. In relation to chemical analysis, for exam- 
ple, an oxygen atom will show the same X-ray absorption in oxygen gas, 
in ice, water, or steam, or in sand. 

The energies of X rays are high enough to make their spectra relatively 
simple. The established ‘‘absorbancy index” for X rays is the mass absorp- 
tion coefficient, tm, which equation 5 defines in measurable terms. This 
coefficient contains the information of interest to the analytical chemist 
who wishes to use X-ray absorptometry in chemical analysis or control. 

Ficure 1, a “log-log” plot, shows how pm varies with wavelength for three 
common metals and illustrates the fact that all X-ray absorbancy curves 
are essentially simple. They consist of large regions in which X-ray ab- 
sorption is continuous, but these regions are separated by sharp discontinu- 
ities (“absorption edges’) where the absorbancy suddenly decreases at a 
“critical” wavelength. These wavelengths are the maximum. capable of 
ejecting a certain kind of electron from the atom. It will be noted that 
the critical wavelengths for both the K and L electrons of lead fall within 
the range of FIGURE 1, but that this range includes only the K wavelength 
for copper. Aluminum shows continuous absorption throughout. Each 
critical wavelength is characteristic of an element, so that the positions of 
absorption discontinuities can be used to identify the elements in a sample 
(qualitative analysis). Generally, however, not all critical wavelengths will 
be equally suited to the identification of a given element. 

The contrast between the discontinuous absorption of radiant energy 


thousands of Angstroms in wavelength, and that of X rays, will: be clear 
from the foregoing paragraph. The contrast between continuous absorp- 
tion in the two cases is no less marked. The mass absorption coefficients 
in regions of continuous absorption can be expressed as the following approx- 
imate function of wavelength \ and of atomic number Z: — : 


Mm = (Z'8)NC/A. (1) 


In this equation, V is Avogadro’s number, A is the atomic weight, and C 
is assumed to be constant in the range between adjacent critical absorption 
wavelengths. The scattering of X rays by electrons gives rise to an additive 
term, usually unimportant except at low Z and X, that has been omitted 
from equation 1. 

Because X-ray absorption is an atomic process, the mass absorption coef- 
ficient of a sample usually stands in a simple additive relationship to those 
of the elements present. If, for example, sample S contains elements A, 


ni and C in the proportions by weight x, y, and (1 — x — y), respectively, 
then 


Mm = Xun + yum + (1 — x — y)uS (2) 


without restriction in relation to the chemical or physical state of S, A, B, 
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or €, For the purposes of chemical analysis or control, the mass absorption 
coefficient of a sample depends only on its ultimate composition, and on 
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Figure 1. Log-log plot showing mass absorption coefficient as a function of wavelength for three common 
metals. Note the characteristic discontinuities separated by the uncharacteristic continuous regions. 
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the wavelength of the X-ray beam, if the beam is monochromatic. With 
polychromatic beams, there may occasionally be departures from equation 2 
(“deviations from additivity”), which are not discussed in this review (bibli- 
ographical reference 6, Figure 4). 
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Types of X-Ray Absorpiometry. X-ray absorptometry may be subdivided 
into absorption spectrometry, absorptometry with filtered beams, and ab- 
sorptometry with polychromatic beams, the order being that of decreasing 
restrictions on wavelength. It will not be necessary to discuss the inter- 
mediate subdivision, for its usefulness can be estimated from that of the 
limiting cases for any particular problem. ; 

In absorption spectrometry, the X-ray beam becomes monochromatic 
before it strikes the sample, or the wavelength distribution of the emergent 
beam is measured by diffraction methods. Because this type of absorp- 
tometry can serve to locate critical absorption wavelengths, it can be used 
for qualitative analysis with the limitations indicated below. The decrease 
in absorbancy at such a critical wavelength will naturally be proportional 
to the amount of the corresponding element in the beam. As will be shown 
later, absorption spectrometry can often establish very precisely this de- 
crease in absorbancy and, consequently, the amount of the element. This 
method is unique among absorptometric X-ray methods in being both 
qualitative and quantitative. Consequently, it is applicable to samples 
about which little or nothing is known. 

When X-ray absorptometry is done without regard to the presence of 
absorption edges, it may usually be considered as a simplified case of ab- 
sorptometry with polychromatic beams, and it is then subject to the funda- 
mental limitation described below. 

Absorptometry with a polychromatic X-ray beam somewhat resembles 
colorimetry with white light. Its most fundamental limitation as an ana- 
lytical tool is implicit in equation 1. Because all atoms absorb X rays, 
such measurements cannot be specific. They cannot ordinarily identify 
unknown elements in a sample, and (as a corollary) they cannot give the 
composition of a sample containing unknown elements, though they can 
usually show whether such a sample has an assumed composition. Figure 1 
abundantly illustrates the lack of specificity inherent in equation 1. Note 
particularly that the mass absorption coefficients for copper and aluminum 
become virtually identical just beyond the K absorption edge of the former. 

Experimental Method and Nature of Results. The measurement of X-ray 
absorption does not differ in principle from similar measurements common 
in analytical laboratories. For example, polychromatic X rays from a 
source (usually a standard X-ray tube) might be passed through the sam- 
ple contained in a cell, whereupon the emergent beam would enter a detector 
to yield an electric current, 7, proportional to the beam intensity, J. When 
absorption spectrometry is being done, there must be provision for suitably 
diffracting the beam before it strikes the sample (to render it monochro- 
matic) or after emergence (to analyze it). Neither operation is possible 
without radically reducing the intensity (say, by several powers of 10) 
below that of the corresponding polychromatic beam. 

Under the simplest conditions, the decreases in beam intensity and in 
output current that measure X-ray absorption by a sample in a cell follow 
the usual exponential law, which may be written: 


log[Zo/I] = logfio/i] = km. (3) 
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In this equation, Zo and iy refer to the empty cell, & is a proportionality 
constant that contains the information of chemical interest, and m is the 
mass of the sample. The effect of certain departures from “simplest con- 
ditions” (“deviations from linearity”) is discussed below. (For other, less 
important departures, see the original literature.*-°) 

By “deviation from linearity” is meant a curvature in a plot of i on a 
logarithmic ordinate against m (or thickness), in other words, a variation in 
k, the slope of such a plot, and a departure from equation 3. Such a devi- 
ation will occur even with a monochromatic beam if the detection and am- 
plification system becomes nonlinear, so that 7 no longer varies directly as I. 
One remedy is to restrict the range of operation until satisfactory linearity 
_isachieved. For example, when the range of sample masses is chosen large 
enough from the viewpoints of convenience and precision, i) may equal 
1000 7. For samples of masses m, and m2, however, equation 2 takes the 
more practical form 


log(Z1/T2] = log[t:/z2] = k[m2 — mi]. (4) 


Another important deviation from linearity peculiar to polychromatic 
beams is a consequence of the fundamental nature of X-ray absorption. 
The proportionality constant, k, is related as follows to the mass absorp- 
tion coefficient of the sample 


Mm = 2.303ak, (5) 


in which a is the cross-sectional area (sq. cm.) of the cell containing the 
sample. For most elements, um is known to increase approximately with 
the cube of the wavelength (equation 1) between adjacent critical absorp- 
tion wavelengths. (The fact that the lines in FicuRE 1 have a slope near 3 
illustrates this point.) Consequently, when a polychromatic beam passes 
through a sample, the longer wavelengths are the more strongly absorbed, 
so that the beam becomes “‘harder”’ as it progresses. It follows that & will 
decrease as m, the mass of sample, increases. This behavior is extremely 
important to the analytical chemist who plans to use polychromatic beams 
for precise work. 

The effect of the deviation just discussed can often be estimated if the 
change in the “effective wavelength” of the polychromatic beam is known. 
(The effective wavelength of a polychromatic beam is defined as that of a 
monochromatic beam absorbed to the same extent under the experimental 
conditions.) Ficure 2 illustrates the points discussed (see bibliographical 
reference 7, pp. 862-3). 


X-Ray Absorption Spectrometry 


It has been mentioned that absorption measurements near a critical wave- 
length can yield both qualitative and quantitative information. Glocker 
and Frohnmayer,’ using photographic methods, did the classical work in 
this field twenty-five years ago with results that were excellent for the times. 

Their method will be illustrated by applying it to determining the pro- 
portion of total bromine (free or combined) x in a sample weighing m grams. 
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Ficure 2. X-ray absorbancy curves for aluminum over two ranges of thickness. Decrease in effective 


wavelength with increasing thickness of sample is principally responsible for curvature shown, which grows 
more pronounced as wavelength increases. 


The experimental data will be the results of absorbancy measurements at 
wavelengths just above and just below (A and 2’, respectively) the critical 
absorption wavelength, 0.918 A. On the basis of equations 2, 3, and 5, 

log Io/I = kpymx + km(1 — x) (A > 0.918 A.) (6) 


log 10/1’ = kama + kgm(1 — x) (A’ < 0.918 A.) (7) 
where 


kar > kpr (see FIGURE 4), (8) 
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In these equations, S refers to the sample less the bromine, and the /,’s 
usually apply to the empty cell. Now, the more nearly identical are \ and 
\’, the more nearly (in the simplest case) 


Ig = Ip and kg = ks. (9) 


If equation 9 is valid, then by substitution, and by subtracting equa- 
tion 7 from equation 6, one obtains 


log I’/I = (kp: — Ra:)mx = —cmx. (10) 


The constant c, which is characteristic for each element and for each 
absorption discontinuity, can be determined empirically or calculated from 
known mass absorption coefficients by equation 5. When either part of 
equation 9 does not hold, the effect is to introduce into equation 10 factors 
that can be determined by making additional measurements. 

Glocker and Frohnmayer determined the characteristic constant c for nine 
elements (bibliographical reference 9, Table 4) ranging in atomic numbers 
from 42 (molybdenum) to 90 (thorium). They proved that identical re- 
sults could be obtained with the sample in the primary (polychromatic) or 
in the diffracted (monochromatic) beam. The method was applied with 
good results to the determination of barium in glass, of antimony in a sili- 
cate, of hafnium in the mineral alvite, and of molybdenum, antimony, 
barium, and lanthanum in a solution of their salts. 

Andrews’ used an ionization chamber as detector to determine iron 
(0.44 per cent) in beryllium according to a procedure resembling that just 
described. Engstrém succeeded in applying the method to microscopic 
sections of tissue and claimed that “in analyses of calcium and phosphorus 
in biological material, quantities of 10~" to 10~'! gram have been determined 
by the method with an error of 10 per cent.”"" These amounts seem very 
small indeed. 

Limitations. The element to be determined may be present in an amount 
so small that the absorbancy of the rest of the sample near the critical wave- 
length is great enough to make the method insensitive (see equations 6, 7, 
and 10), or the characteristic constant c may not be large enough to permit 
attainment of the precision desired. Measurements at certain critical ab- 
sorption wavelengths, e.g., the K series for the light elements, are difficult 
to make, because X rays of long wavelength are strongly absorbed by air 
and by window materials. There may be interference owing to the near- 
ness of critical wavelengths characteristic of other elements in the sample, 
e.g., the K wavelengths of lead and thallium differ by only 0.004 A. The 
seriousness of these limitations can sometimes be reduced by selecting a 
more favorable critical absorption wavelength. More important, however, 
is the promise of progress in this direction that is held out by recent im- 
provements in the equipment available for such work. 

Dow Spectrometer. An outstanding example of such equipment is the 
Dow automatic X-ray absorption spectrometer being developed by Frevel 
and North,” who are now preparing to publish an account of their work. 
The front of the instrument is shown in FIGURE 3. The spectrometer itself 
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is housed in the case that supports control and recording apparatus and 


extends a considerable distance beyond the right edge of the photograph. 


‘ 


Ficure 3. Dow automatic X-ray recording spectrometer, front view. Note cell on spectrometer case, 
which extends off photograph to right. Driving mechanism (left front) and vertical X-ray tube (center) 
are shown. (Courtesy Dow Chemical Company.) 


A cone of polychromatic X rays (generated by the tube mounted vertically 
in front of the spectrometer case) passes through the cell containing the 
sample and into the case, where the unabsorbed radiant energy strikes a 
multiple-crystal “lens,” comprised of four sodium chloride crystals, whose 
monochromatic beams are focused on a Geiger counter so that the sum of 
their intensities can be automatically recorded as an output current. Vari- 
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ation of wavelength is accomplished by a lathe lead screw (see the front of 
FIGURE 3) which moves the lens and the detector (at twice the lens speed) 
along the optical axis. “Voltage is stable to within 0.05 per cent during 
measurement.... The probable error of the mean of a single intensity 
measurement is 0.5 per cent for all measurements made for a period of 
approximately five minutes at rates of 10,000 counts per minute and higher. 
Repeated measurements give agreement within one per cent. The time 
for an I/Ig measurement is 10 to 15 minutes.’”’” 

It seems safe to predict from these data that the Dow spectrometer will 
make possible greater accuracy than was achieved by Glocker and Frohn- 
mayer, who found, for example, 5.45 per cent barium in a glass that gave 
5.8 per cent upon chemical analysis. The increased convenience in meas- 
urement arising out of the use of modern instrumentation, however, will 
often outweigh any such gain in accuracy. In the barium determination 
just mentioned, Glocker and Frohnmayer used a 90-minute exposure, after 
which the film had to be painstakingly developed and microphotometered. 
On the Dow instrument, output currents proportional to intensities are 
simply taken from a chart (bibliographical reference 5, Figure 3). 

Ficure 4 shows how experimental! data for the determination of bromine 
(cf. the first part of this section) on the Dow spectrometer are extrapolated 
to give the decrease in absorbancy at the critical wavelength. 

As regards precision, time required, and range of determinable elements 
(now down to atomic number 22), Frevel and North have improved the 
method of Glocker and Frohnmayer to the point where it deserves serious 
consideration by analytical chemists. Analyses for the lighter elements can 
of course be carried out on the Dow instrument by making measurements 
at various wavelengths in regions of continuous absorption. 


X-Ray Absorptometry with Polychromatic Beams 


The use of monochromatic X-ray beams is desirable in chemical analysis 
because it simplifies the interpretation of results, but there are many appli- 
cations in which the concomitant reduction in beam intensity cannot (or 
need not) be tolerated. Outstanding among these are control operations, 
two of which will be described below. 

Absorptometry with polychromatic beams received a powerful impetus 
from the invention of the multiplier phototube!: !° and its subsequent con- 
version into a photoelectric X-ray detector by the addition of a suitable 
phosphor (see FIGURE 5). Morgan,'® who had been investigating photo- 
electric cells for the control of roentgenographic exposures in medical work, 
was apparently the first to make this conversion. Independently and 
somewhat later, Smith” and Moriarty constructed such a detector in 
solving a war problem that was fundamentally identical with many prob- 
lems in chemical control. 

Control Problems. If a hand-grenade fuse contains too little powder, the 
grenade explodes prematurely and may endanger the man who is trying to 
use it. An infallible, non-destructive, rapid test is obviously required, and 
none was available until Smith and Moriarty provided the one described 
below. 
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ere mounted upright on a belt that carried them 
beam. The absorbance of a fuse contain- 
ing the proper amount of powder was great enough so that the intensity of 
the transmitted beam, as measured by a photoelectric X-ray detector, was 
too low to trip the detector circuit. With a defective fuse in the beam, 
however, the increased output current from the multiplier phototube set 
into operation four means of identifying the defective fuse. The equipment 


Fuses to be tested w 
through a polychromatic X-ray 


-3 -.2 -.| @) | 
LOGe (A) 


Ficure 4. Extrapolation of experimental data in region near absorption edge to give bromine content 
of sample. A log-log plot permits linear extrapolation (cf. FIGURE 1). (Courtesy Dow Chemical Company.) 


inspected fuses at a rate of approximately 4,000 an hour, and this rate was 
not limited by the detecting or recording apparatus. 

Clapp and Pohl'* solved another important control problem by using a 
polychromatic X-ray beam to measure the thickness of steel strip. During 
the measurement, the strip is hot (1400° to 1700°F.), moving (2000 feet per 
minute horizontally with possible vertical vibrations up to several inches in 
amplitude), and subjected to a spray of cooling water. As is being done to 
an increasing extent in instrumental analysis, the measurement is accom- 
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plished by means of a servo system, which is adequately described in bibli- 
ographical reference 19. 

It seems reasonable to hope that X-ray absorptometry with polychro- 
matic beams can be used to solve various other production problems, such 
as assessing roughly the quality of crushed minerals on a conveyer belt, or 
controlling the addition of a material, e.g., tetraethyl lead to gasoline, to a 
moving ‘‘base stock.” 

Chemical Analysis. The use of polychromatic X-ray beams in chemical 
analysis is not new. Fuller”? used photographic means to compare the 
absorbancies of alloys for the purpose of establishing their composition. 
Aborn and Brown,” acting upon a suggestion of Calingaert, applied absorp- 
tometry of this type to the determination of tetraethy] lead in gasoline and 
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Ficure 5. Schematic diagram of photoelectric X-ray detector. Phosphor converts X-rays into light 
that liberates electrons from photocathode O. These are guided electrically to the successive dynodes 1 to 
9, where multiplication to give secondary electrons occurs. The greatly (10° to 106 times) amplified beam is 


finally gathered by anode 10 for subsequent external amplification (over 104 times in photometer of FIGURE 
6), if desired. 


used an ionization chamber as detector. Recently, Sullivan and Friedman” 
employed a Geiger counter for the same purpose. 

The successful solution of the fuse-testing problem led to an investigation 
of the photoelectric X-ray detector as a tool for chemical analysis.” The 
laboratory photometer incorporating this detector is shown in FIGURE 6. 
Early work on solids, liquids, and gases was done by the direct method.’ 
The intensity of the polychromatic X-ray beam was adjusted to a standard 
initial value by varying the X-ray tube voltage until the desired output 
current was obtained with a standard thickness of aluminum in the beam, 
the voltage across the detector and the amplifier setting being fixed. Out- 
put currents obtained with known weights of sample in the beam were then 
used to give information about the composition of the sample. In general, 
these output currents are plotted as in FIGURE 2, and the slopes of the 
curves (or straight lines, equation 4) are interpreted according to equations 


5 and 2. 
To illustrate such an interpretation, a common analytical problem, the 
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determination of chlorine in a chlorinated hydrocarbon polymer, a be 
briefly discussed.’? In such a polymer, chlorine may be thought of as re- 
placing hydrogen. Because the former element has by far the a. mass 
absorption coefficient, this replacement will increase k, and this increase 


Ficure 6. Laboratory X-ray photometer. A. Photoelectric X-ray detector, 
ple. D. CA-5 X-ray tube and housing. E. Milliammeter. 
Regulated power supply for amplifier tubes and phototube anode voltage. 


B. Sample cell. C. Sam- 
F. Amplifier and rectifier vacuum tubes. G. 


H. Control panel. 


will be proportional to the extent of the replacement. Consequently, a 
plot of k against chlorine content for a family of such polymers should 
show an intercept on the & axis that represents k for the unchlorinated 
hydrocarbon. The other k values should fall upon a straight line whose 
slope depends upon the difference between the mass absorption coefficient 


Liebhafsky: X-Ray Absorption 1009 


for chlorine and that for hydrogen. FicurE 7 shows that these expectations 
are realized. 

The abscissas in FIGURE 7 are the results of chlorine determinations by a 
conventional method. The solid line shows that X-ray and conventional 
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Ficure 7. Chlorine contents of chlorinated hydrocarbon polymers. A comparison of results by the 
direct method with those from conventional analyses. 


results are in good agreement. It is obvious that these data could be used 
as an empirical calibration curve, with the help of which chlorine deter- 
minations could be carried out on unknowns of this type about as satis- 
factorily as by the conventional method, and in one tenth the time. The 
circles and the dotted line in FicuRE 7 are calculated results based upon 
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known mass absorption coefficients (bibliographical reference 7). The slope- 
intercept relationship of FIGURE 7 is common to many other problems, such 
as the determination of sulfur or of additives (such as tetraethyl lead) in 
petroleum products. 
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y the comparative method (vari- 

The direct method has been used in this laboratory in point-to-point 
explorations of impregnated materials to test their uniformity. For this 
work, the diameter of the X-ray beam was reduced (in some cases to 0.15 
cm.), and considerably over ten thousand output-current readings were 
taken. Each reading gave information in a matter of seconds that would 
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have required several hours to obtain by conventional methods. Had con- 
ventional methods been used, moreover, analytical and performance data 
could not have been obtained on the same specimen, an overriding consid- 
eration in this case. 
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FicureE 9. General Electric X-ray photometer. Front view, oblique from left, with parts named. 


Though the direct method is well suited to exploratory work and yields 
results precise enough for many purposes, it is subject to uncertainties 
arising from voltage fluctuations and from changes in the effective wave- 
length of the polychromatic beam.® For an unfiltered beam of this kind, 
the variation of output current with primary voltage is much more pro- 
nounced than for a monochromatic beam, because the change in the total 
X-ray output affects the current readings in the former case. Rapid com- 
mutation in the beam between the unknown and a suitable standard can 
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Ficure 10. Control charts to illustrate voltage fluctuation problem. Ten readings constitute each 
group. I. Output currents on laboratory photometer (r1cuRE 6). II. Thicknesses of aluminum required 
to balance 250 mils of aluminum against 280 mils (both approximate) on the G. E. X-ray photometer. Regu- 
lated voltage source. III. Experiments of Chart II with 110-volt line as voltage source. 
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reduce the uncertainties due to both causes, and the comparative method 
of X-ray absorptometry uses this procedure. Readings of output currents 
are taken alternately for standard and unknown. From these readings, 
the amount of standard equivalent in absorbancy to the unknown can be 
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calculated. As a consequence, the interpretation of the results does not 
directly involve output currents, and this is the great advantage of the com- 
parative method. In many cases, aluminum is a satisfactory standard. 

The laboratory photometer was used in carrying out the following three 
types of determinations by the comparative method’: identifications of cer- 
tain new compounds, determination of tetraethyl lead fluid in gasoline, and 
determination of sulfur in crude oil. The results were generally satisfac- 
tory, and the expected improvement over the direct method was realized. 
Frcure 8 shows a set of such results for one of a number of similar problems 
in the petroleum industry. The measurements were made on unknown 
samples kindly sent us by the Ethyl Corporation; details of the work have 
recently appeared.® 

General Electric X-Ray Photometer. When the direct method is being 
used with polychromatic beams, an extremely constant voltage source is 
necessary if the results are to be precise. This point is illustrated for the 
laboratory photometer (FIGURE 6) by the upper control chart™ in FIGURE 10. 
Although the voltage source was a synchronous motor-generator set, the 
output of which was further regulated by two Sola transformers in series, 
it is clear that the output current readings in this chart show an “‘assigna- 
ble cause of variation,” which must have been the inconstancy of the voltage 
source. The introduction of the split-beam technique and the comparative 
method should improve this situation. Michel and Rich,” in developing 
the G. E. X-ray photometer, decided to use a “‘chopper’’ on a single beam 
in order to accomplish this purpose, but other methods are conceivable. 

In the G. E. X-ray photometer (FIGURE 9), the rate of synchronous com- 
mutation between standard and sample is so high (30 times per second) 
that the ordinary 110-volt alternating current usually requires no regula- 
tion to be satisfactory as a power source. ‘This point is illustrated by the 
two lower control charts in FIGURE 10. The regulated voltage source of 
the middle chart was that used for the laboratory photometer. The two 
lower charts show systems much more nearly “in control” than that of the 
upper chart, and they indicate that our 110-volt line is a satisfactory power 
source for the General Electric X-ray photometer. 

The dotted lines in FIGURE 10 are the conventional 3 ¢ limits (biblio- 
graphical reference 24, p. 53, §[ 14) from which the precision may be esti- 
mated when the system is in control. In the middle chart, for example, 
the total thickness of aluminum in the beam (ca. 280 mils) should be estab- 
lished to within approximately 0.3 mil by the average of virtually each 
group of ten readings. 


Conclusion 


The Geiger counters and photoelectric X-ray detectors now available have 
increased the usefulness of X-ray absorptometry. To a large extent, these 
devices are complementary. At present, Geiger counters are favored in 
X-ray spectrometry, although the photoelectric detector may eventually 
enter this field.2 Absorptometry with polychromatic beams, especially in 
control applications and point-to-point explorations, where speed is vital, 1s 
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best done with the photoelectric detector, because it is virtually in- 


stantaneous. : ; ; é 
Geiger counters are incorporated in the improved diffraction equipment 


(Norelco, G. E. XRD-3) now available, which makes this equipment readily 
adaptable to X-ray absorptometry. In fact, it would be logical for absorp- 
tion methods to be tested in research laboratories on such equipment before 
they are applied in the plant on commercial instruments such as the G. E. 
X-ray photometer. Recently, Kehl and Hart” have determined sulfur in 
hydrocarbons by X-ray absorptometry on just such a diffraction apparatus, 
and Vollmar, Petterson, and Petruzzelli® presented work of the second kind 
at the San Francisco meeting of the American Chemical Society. 
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THE ANALYTICAL SPECIFICITY OF COUNTERCURRENT 
DISTRIBUTION 


By J. Delafield Gregory and Lyman C. Craig 
The Rockefeller Institute for Medical Research, New York, N.Y. 


The rapid development of modern biochemistry has brought the realiza- 
tion that classical techniques of isolation and purification are frequently 
inadequate to cope with the types of compounds encountered. Instability, 
high molecular weight, and annoying physical properties, such as failure to 
crystallize, are often handicaps, together with the likelihood that, in a bio- 
logical environment, closely related analogs and isomers may be found 
occurring together to aggravate the problems of separation. The need for 
analytical methods, which can be applied as generally as possible in such 
cases, has led to the development of countercurrent distribution in this 
laboratory. Since it is essential to have techniques which are sensitive and 
specific, but which operate under the mildest experimental conditions, the 
method has been designed to utilize the advantages of working with dilute 
solutions. 

- In this procedure, a compound or mixture is partitioned between series 
of immiscible phases in a simple countercurrent order in separate contact- 
ing cells, every individual partition attaining essentially complete equi- 
librium. Only the case of partition between two liquid or solvent phases 
will be considered here, although the principles apply regardless of the na- 
ture of the phases.! 

Before considering the potentialities of the method as an analytical aor 
- it will be necessary to outline briefly the actual procedure and some of the 
simple mathematics which furnish its basis and validity. Then, the appa- 
ratus currently in use and the type of results which can be obtained will be 
described. 

When a single pure solute is dissolved in either one of a pair of immiscible 
phases and the phases are shaken together, the solute will partition itself 
between the two in a characteristic proportion at equilibrium.” The ratio 
of the fraction in the upper phase to that in the lower (assuming equal 
volumes) is called the partition ratio or distribution coefficient, K. The 
following development deals only with sufficiently dilute solutions, so that 
the partition ratio is not significantly affected by the total concentration, 
and so that, in case of a mixed solute, every individual species behaves 
independently of the others and has its own reproducible K. This ratio, 
it should be emphasized, is not directly related to the relative solubilities 
in the two phases, since solubility implies saturated solutions where condi- 
tions are very different from those in the dilute solutions considered here. 

The case of a pure solute which partitions equally between the two phases 
is taken for illustration in FIGURE 1, assuming equal volumes of the phases 
throughout. After the first equilibration, each phase contains half the 
material. They are separated and the upper transferred to the adjacent 
cell, which contains fresh lower phase; fresh upper is added to cell 0. After 
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equilibration of both cells, the uppers are each moved one cell to ay right 
and the equilibration repeated. Every time the whole series of upper 
phases is moved one stage, one transfer has been made. At any given 
stage, the numbers in FIGURE 1 represent the fraction of the original start- 
ing solute now found in each phase. The overall effect of progressive 
transfers is to cause a series of upper phase portions to pass over a series 
of lower phases. For the sake of consistency, it will be considered that 
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Ficure 1. Schematic diagram of countercurrent distribution of a substance with partition ratio of unity. 


the upper phases always move to the right. The two series “flow,” although 
discontinuously, countercurrent to each other. 

When a similar analysis is carried through, using a general expression for 
the fractions in the upper and lower phases, instead of this specific case, it 
is found that for any number of transfers, ”, the contents of each of the 


n+ 1 cells at the end of the distribution is given by the terms of the ex- 
pansion of the binomial? 
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For a given number of transfers, when the terms of the expansion are plotted 
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as ordinates against their position in the series of cells, a curve such as one 
of those in FIGURE 2 results. The shape and position of the curve depend 
on the partition ratio of the solute, since the upper phase, moving to the 
right, carries it along further the higher its K, that is, the more it favors 
the upper phase. The curve of a substance with K greater than unity is 
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Ficure 2, Theoretical curves for twenty-four transfer distribution. 


a mirror image of the curve for the reciprocal K, but at the other end of 
the series. 

When the starting material contains more than one species, each will 
form its own particular distribution, and the curve found at the end of the 
process by the determination of the material in each tube will be the sum 
of the individual binomial curves. 

If a particular K is considered (equal to unity for illustration), and the 
other variable, the number of applied transfers, is changed, the effect can 
be seen in FIGURE 3. With more transfers, the curve becomes lower and 
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broader, but its width does not increase as fast as the total number of tubes, 
as shown in the table in FIGURE 3. The ability to apply larger numbers of 
transfers allows more complex mixtures to be fractionated, since there is 
room in the distribution pattern for more discrete curves to appear. 

For a given K, the curve always maintains the same position relative to 
the total number of tubes (the center when K = 1). This follows from the 
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Ficure 3. Theoretical curves for K = 1 and various numbers of transfers, 


relationship of V, the number of the tube containing the maximum, to a 
the number of transfers‘ (the total number of tubes is + 1). 


N= n( es), (2) 


apes continuous function is an approximation as applied to the discontinu- 
ous distribution process and is not satisfactory for » less than about twenty 
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or for values of K far from unity. Under those conditions, more exact and 
cumbersome calculations based on the binomial expansion must be used. 
The position of the peak of a pure substance in a distribution is the most 
accurate measure of its partition ratio, much more so than any determina- 
tion based on a single partition. 

The available number of contacting cells limits the number of transfers. 
which can be done in this simple way. But the operation can be carried 
further by an extension of the procedure, which is called ‘withdrawal.’ © 
After the cells are filled, following the scheme of FicurRE 1, further transfers 
give upper phases emerging at the right hand end of the series, while fresh 
upper phases are fed in from the left at each stage. The effect produced is 
that of an unlimited series of upper phase portions passing over a fixed and 
limited set of lower phases. Similarly, a series of lower phases can be made 
to pass under the fixed uppers, or the two effects can be combined to allow 
both phases to be withdrawn on alternate transfers. Since each withdrawn 
fraction has experienced a different and progressively increasing number of 
transfers, the curve of a pure solute, when it passes into the withdrawn 
series, is not quite symmetrical. The situation can be dealt with mathe- 
matically by simple approximations, and the withdrawal curves can be 
interpreted as readily as the curves in the basic part of the distribution. 
This aspect will be treated fully elsewhere. 

For distributions of twenty-four transfers or less, and, in any case, where 
part of the solute remains in the first or last tube, so that its curve is not 
symmetrical, the exact formula of the binomial expansion must be used in 
constructing curves. The fraction of the solute in any tube 7 can be found 
from formula 3, and there are many other relationships which simplify the 
calculations.’ 
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On the other hand, for large numbers of transfers, the arithmetic becomes 
unnecessarily cumbersome. The exact binomial curves then approach nor- 
mal distribution curves so closely that the mathematics of probability can 
be substituted as a very excellent approximation. The only formulas 
needed for calculating a theoretical curve are the following: 
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The fraction of the solute in the maximum tube is given by yo, and yz rep- 
resents the fraction in a tube, x tubes distant from the maximum on either 
side. Only enough values of yz, suitably spaced, need be calculated so 
that the curve can be drawn readily. 

Another useful quantity, stemming from the mathematics of probability, 
is the standard deviation, c. A range of + 3 o about the maximum in- 
cludes more than ninety-nine per cent of the area of a curve,® and, for 
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practical purposes, these are considered the limits of the distribution. In 
terms of the symbols used in these partition studies, 
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For many preliminary calculations and estimations of the outcome of dis- 
tributions, the widths and positions of the curves as obtained from equa- 
tions 2 and 5 are sufficient. 

Even though the experimental and calculated points are connected by 
smooth curves, and continuous functions are used for approximation, it 
should be emphasized that the operation itself is discontinuous, the only 
real values on the graphs being the points themselves. 

Calculated curves have two chief uses, one of which is the construction 
of hypothetical situations which serve as guides in selecting conditions for 
actual experiments. In addition, it is possible to interpret the results of 
distributions by fitting calculated curves to them and determining the 
minimum number of components required to account for the experimental 
patterns. A distribution of a mixture of four fatty acids’ is shown in 
FIGURE 4. The resulting curve shows four peaks, which can be closely fitted 
by the sum of the four theoretical curves plotted on the same graph. ' The 
validity of the analysis in this case was confirmed by isolation and redistribu- 
tion of the material in tubes 15, 16, and 17. A curve resulted which was 
accurately matched by a single theoretical curve. The quantitative aspects 
of this curve-fitting procedure, when used as a method of determining 


purity, will be discussed later. “i 

The application of these principles as an analytical method is. justified as 
follows: A single partition between two phases is sufficient in theory to cause 
a partial separation of a mixed solute if there is any difference whatever in 
the partition ratios of the components. Simply determining the mass of 
material in each phase is not revealing at this stage. It is necessary to use 
a more specific analysis which can, in its own right, differentiate between 
thé components of the mixture. If more partitions are done; ‘that’is, if a 
distribution is made, then the mere determination of weight in the' various 
fractions is enough to establish whether species with different’ partition 
ratios are present. In this way, weight, the most fundamental and im- 
partial measurable property, is made to give very specific analytical infor- 
mation. Weight is the preferred analytical unit for interpreting most dis- 
tributions, especially for exploring new situations where some component 
of a mixture might be overlooked if a more specific analytical method such 
as spectrophotometry, titration, or bioassay were relied upon. 

In this or any other partition method, the sensitivity and separating 
power are functions of two experimental conditions: the number of effective 
or equivalent partitions which can be applied, a mechanical problem, and 
the obtainable difference between the partition ratios in a mixture. The 
latter involves the careful study of solvents and their effects. The greater 
the number of transfers in a distribution, the more individual operations 
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contribute to the final result, which then represents the average of a larger 
number of single partitions. This minimizes the effect of incidental errors 
in transfers and manipulation. For example, in a routine distribution using 
fifty-four cells and four hundred transfers, well over 20,000 useful extractions 
are performed. It is evident that the experimental curves bear a real rela- 
tionship to the normal distribution or probability curves on which the 
mathematical treatment rests. 

Consequently, one aspect of the subject which has received much atten- 
tion in this laboratory is the design of apparatus to simplify making large 
numbers of transfers. Stainless steel countercurrent distribution instru- 
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Ficure 4. Twenty-four transfer distribution of a mixture of (left to right) acetic, propionic, butyric, 
and valeric acids. ©, experimental curve; @, theoretical curves. 


ments have been in use for some time.‘>* The one shown in FIGURE 5 
contains fifty-four tubes arranged in a circle. The open tops and bottoms 
are all closed simultaneously by two large glass discs, and a sliding joint 
between the upper and lower sections of the tubes allows all the upper sec- 
tions to be shifted as a unit with respect to the lower ones about the central 
shaft. The lower phase exactly fills each lower section, and the upper 
phases are carried around from one to the next. At each stage, the solvents 
are thoroughly mixed by tumbling the whole assembly on the side bearings, 
the equivalent of shaking fifty-four separatory funnels all at once. 
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With a fifty-four tube instrument, up to three or four hundred transfers 
can be made profitably using one of the withdrawal procedures. Such a 
run, together with the analysis of the results by evaporating aliquots and 
weighing the residues, requires three to four man-days of work. There are 
many fractionation problems for which this would be more than ample, but 
as more difficult situations are undertaken and more separating power is 


Ficure 5. Fifty-four tube stainless steel distribution instrument. 


needed, methods are being developed for making still more transfers quickly 
and easily. An all-glass apparatus, operating on a simple principle of de- 
cantation, has been designed® and has already established its usefulness.* 
A single unit or tube is seen in FIGURE 6. The equilibration is made in 
chamber A by rocking the tube about 30° above and below the horizontal. 


* A more recent design may be found in: Craic, L. C. 1950. Anal. Chem. 22: 1346 
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After the phases have separated, the unit is tilted backward 90° so that all 
the upper phase runs out into the auxiliary transfer chamber, C. The 
lower phase remains in the end section of A, which it fills exactly. When 
the tube is righted, the upper phase runs out into the next unit in line. 
Extension of the train is virtually unlimited, as any number of units may be 
set up side by side. 

FicuRE 7 shows an assembly containing one hundred and eight tubes. 
Transfers can usually be made at the rate of thirty to forty per hour, which 
means about 4,000 individual extractions per hour. The advantages of 
glass construction and the absence of moving or sliding parts are obvious. 
The assistance of Mr. Otto Post has been invaluable in the design and 
construction of these instruments. 


FicurE 6. Units from glass distribution instrument. 


It is interesting to compare certain aspects of such an apparatus with a 
continuous liquid-liquid extraction column. With regard to rates, when 
solvent is being fed in at one end of the train and extract withdrawn from 
the other, the throughput may be more than 400 ml. per hour. There are 
two important respects, however, in which the method differs from any 
continuous process: this distribution operates unquestionably at equilibrium, 
and the exact number of stages (which are practically theoretical stages) is 
known, so that no assumptions are necessary for the mathematical inter- 
pretation. 

The practical working range of partition ratios depends on the purpose 
of the distribution and on the number of transfers which can be applied to 


1024 Annals New York Academy of Sciences 


to the solute. The maximum separation per transfer of a binary ser 
occurs when the partition ratios of the two species are reciprocals of eac 

other.?: 8 Even though the curves are broadest when the K S are ae 
unity, transfers have the greatest effect on the solutes in this region, oe 
differences in their properties are exploited to the utmost. For genera 
preparative purposes, however, in a small distribution of twenty-four trans- 
fers, the working range of K might be approximately 0.2 to 5, for one 
hundred transfers 0.1 to 10, and for one thousand transfers 0.01 to 100. 
Although efficiency is lost when the K is either high or low, the fact that 
many transfers can now be made quickly offsets the disadvantage to some 
extent and makes it practical to handle a much wider variety of solutes for 
which systems have not been found giving K’s in the optimum range. 
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FIGURE 7. 108-Tube glass distribution instrument. 


Many amino acids and peptides, for example, can be distributed between 
an aqueous phase and a partially miscible alcohol,!! even though most of 
them greatly favor the aqueous phase. 

Several questions arise in connection with the applicability of counter- 
current distribution. How many transfers are necessary to separate a given 
pair of solutes to a given degree of completeness? What is the smallest 
proportion of a given impurity which can be detected under conditions 
likely to be encountered? What is the relationship between the nature of 
a solvent system and its ability to separate particular types of mixtures? 
The mathematics of the binomial expansion answers the former points. 
The problem of the specificity of solvents is not so simple, and some gen- 
eralities and intuitive ideas must suffice for the present. 

In order. to represent the separations which can be attained, calculated 
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curves have been set up for various situations. For the sake of simplicity, 
a mixture of two species in equal amount may be assumed, distributed under 
the most favorable conditions. This means that the partition ratios are 
reciprocals, the resulting curve being symmetrical about the center of the 
whole distribution. If a degree of separation is arbitrarily chosen such that 
fifty per cent of each component is recovered essentially pure, the distance 
between the peaks is equal to one half the width of either curve, or 3 o, as 
in FIGURE 8. The separability or “separation factor” of the two compo- 
nents may conveniently be represented by 8, the ratio of the larger K to 
the smaller. In this special case, 8 = 1/K*. This quantity is analogous to 
a, the term used in distillation theory for the ratio of the vapor pressures 
in a binary mixture. When 6 = 1, separation is impossible in the given 
solvent system, and, as #8 increases, separation becomes easier. The rela- 
tionship between 8 and the number of transfers required to achieve this 
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Ficure 8. Relationship between separability of a binary mixture and number of transfers required for 
arbitrary degree of separation. 


particular spread may be derived as follows (the positions of the maxima 
are given by formula 2): 


K ae NT ( 1 ) 
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The distance between the maxima is: 
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9K 
whence n(1 — K) = 3-\/nK and n= G- e (8) 
From this relationship and the graph in FIGURE 8, it is seen that the 
required number of transfers increases very rapidly as 6 decreases in the 
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region between two and unity. No matter what the 8, if the postulated 
degree of separation is attained and all the material is recovered in two 
fractions by dividing at the minimum of the curve, each part will be better 
than ninety-three per cent pure. If essentially complete separation were 
demanded, the distance between the peaks would have to be 6 a, with cor- 
responding changes in the formulas.” 

When a distribution is used for preparative isolation, it is usually desira- 
ble to start with as much solute as possible. The danger then arises of 
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FicureE 9. 100-Transfer distribution; effect of scattering solute initially in seven tubes 


overloading the solvent system to the point where the partition ratio is no 
longer sufficiently constant and the resulting curve is distorted or skewed 
This may be largely offset by scattering the solute initially in several ad- 
jacent tubes. Ina one hundred-transfer distribution, if the solute is started 
in seven tubes, the result will be the sum of the seven small overlapping 
curves in FIGURE 9, given by the solid line. If the same amount of solute 
were started all in one tube, the dashed curve would result. Unless the 
ultimate in accuracy is required, the difference in these results may be dis- 
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regarded, and the advantages of beginning with only one seventh the con- 
centration may be very great. The discrepancy will be no worse than this 
if the material is scattered in a number of tubes equal to about twenty per 
cent of the ultimate width of the curve at the end of the distribution. When 
K = 1, one hundred transfers gives a curve thirty tubes wide, which would 
allow starting in six or seven tubes. 

So far, this discussion has been from the point of view of separation and 
isolation from gross mixtures. The study of purity is essentially the same 
problem, with the understanding that contaminants are usually in smaller 
proportion and that emphasis is on the analytical results. Since the abso- 
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lute purity of a substance can never be proved and can be confirmed only 
by continued failure to find impurity, the definition of purity must be an 
operational one in terms of the specific techniques involved. Although a 
great deal may be accomplished on a preparative scale, even when the par- 
tition isotherm is not linear and the curves are slightly skewed, for sensitive 
detection of small amounts of impurity it is essential to work at low concen- 
trations and under the most ideal conditions attainable. No distribution 
can reveal a mixture if the components have identical K’s, or if the amount 
of impurity is below the limit of sensitivity of the measurements used to 
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obtain the distribution curve. Identical partition ratios are certainly un- 
likely (except for optical enantiomorphs) especially if more than one solvent 
system can be applied, but it is important to know exactly what degree of 
sensitivity may be expected. 

In order to form an idea of the limits to which the technique can be 
carried, the situation in FIGURE 10 may be considered. An impurity of 
ten per cent should be detectable, using one hundred transfers, if the sol- 
vent system is selective enough to give a 6 of 1.2. The outer curve is the 
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Ficure 11. Hypothetical 1000-transfer distribution of a mixture of 90 per cent compound A (K = 1.0) 
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8, per cent compound B (K 1.1). C =A-+ B. D = theoretical curve for single substance fitted 


one which would be found on weight analysis of such a distribution. With 
only this information, a theoretical curve would be calculated to fit as closely 
as possible. The experimental curve deviates most from the theoretical in 
the region of the impurity, and the difference here is significant provided 
the weight determination has a sensitivity at least as good as one per cent 
of the maximum weight. 

If the 8 is lowered to 1.1, however, as in FIGURE 11, the number of trans- 
fers needed to show an impurity of ten per cent approaches one thousand 
In general, one might expect to find such a low 6 only in the case of very 
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closely related isomers, or else as a coincidence which might be eliminated 
by changing to a different system. If the B- values were higher in these 
examples, of course, impurities present in much smaller amount would be 
observable. A valuable characteristic of countercurrent distribution in 
purity studies is that the impurities receive the same treatment as the 
major component and are often as readily isolated.¥ 

The sensitivity of the method as outlined here is based on analysis by 
weight alone. This quantity should always be used when possible because 
it is all-inclusive, but other methods which differentiate between the com- 
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FicureE 12. Isolation of a penicillin from a crude extract of a culture medium. 


ponents can supplement weight and greatly increase the sensitivity of the 
results. A specific analysis such as bioassay is often useful in the isolation 
of natural products, as shown in FIGURE 12, which represents the separation 
in a single operation of a penicillin from a crude extract of the culture." 
With respect to solvents, the only basic requirement is a pair of phases 
between which the material to be investigated has a reasonably constant 
partition ratio in the working range at a suitable temperature. The only 
restriction in the nature or number of solvents which make up the system 
is that two phases must result, when they have all been shaken to equi- 
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librium, which do not cause transformation of the solute and which do not 
interfere with the determination of the results, or, in preparative work, 
with the isolation of the fractionated products. The whole procedure can 
be done in the cold room if it is necessary to protect labile substances, or 
at any other particular temperature which may have been found to give 
useful solvent properties. 

In general, it may be said that the greatest specificity is attained when 
the two phases show the strongest competition for the solutes. Also, the 
phases should be as different as possible, perhaps as measured by their rela- 
tive dielectric constants. Specificity is generally, though not always, re- 
duced when solvent combinations which are nearly miscible are used. In 
the latter case, the composition of the phases is more susceptible to distor- 
tion by the presence of the solute and by small temperature changes. In 
dealing with acidic or basic solutes, a system composed of an aqueous buffer 
and an organic solvent is often very effective,® since it takes advantage of 
subtle differences in acidity or basicity, and since the partition behavior 
can be regulated conveniently through the pH of the buffer. 

In this brief discussion of countercurrent distribution, the analytical speci- 
ficity of the method has been shown to depend on the following factors: 
(a) the number of transfers which can be practically applied (at the present 
time, apparatus is being developed which greatly increases this number); 
(b) the selectivity of the solvent systems, a field in which empirical data 
and experience are the best guides so far; (c) the ideality of the solvent 
systems, so that partition ratios remain constant, giving results which are 
predictable by simple mathematics and without assumptions; and (d) the 
sensitivity and specificity of the measurements used to evaluate the results 
of the distributions. 
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METHOD AND APPLICATION OF DIFFERENTIAL 
THERMAL ANALYSIS* 


By Ralph E. Grim 


Division of Clay Resources and Clay Mineral Technology, State Geological Survey, 
Urbana, Illinois 


The method of differential thermal analysis determines, by suitable appa- 
ratus, the temperature at which thermal reactions take place in a material 
when it is heated continuously to an elevated temperature and, also, the 
intensity and general character of such reactions. Because of dehydration, 
oxidation reactions, destruction of lattice structure, change in crystalline 
phase, efc., many substances undergo thermal reactions when they are 
heated. The temperature of the thermal reactions, their intensity, and 
their general character are characteristics of a particular material and there- 
fore provide criteria for its identification. Further, since such thermal reac- 
tions result from changes within the material heated or from other reactions 
between components of a mixture, the differential thermal method affords 
a means of studying high-temperature phase changes or reactions. 

Thermal reactions of very low intensity can be utilized in the method if 
they begin abruptly and are completed in a short temperature interval. 
The results are recorded generally in the form of a continuous curve on 
which endothermic reactions are conventionally deflected downward and 
exothermic reactions upward from a horizontal base line. 


History of Method 


In an investigation of the components of clays in 1887, Le Chatelier’ 
appears to have set up, for the first time, simple apparatus to measure and 
record thermal reactions taking place in a material when heated and to 
have used these reactions as analytical criteria for studying the components 
of a material. Following Le Chatelier, thermal reaction studies of clays 
were made by several investigators.!: *!, 3-44.45 Jn these early studies, the 
procedure generally used was to place the material in a small platinum 
crucible with a single thermocouple junction in the center of the material. 
The whole mass was placed in a furnace and heated at a rapid and rela- 
tively uniform rate. The thermocouple was attached to a galvanometer 
which was read or recorded photographically at short intervals. Thermal 
reactions in the material caused variations in the galvanometer record as 
compared to that obtained when the furnace contained no sample. The 
record obtained showed the thermal reactions of the material superimposed 
on the heating curve for the furnace alone (FIGURES 1 and 2). In some of 
this early work,*® the method consisted essentially in taking the time-tem- 
perature record of a thermocouple embedded in the specimen and compar- 
ing it with the record given by a second couple placed in a reference sub- 
stance alongside it. In still other*” early work, the so-called Saladin-Le 
Chatelier apparatus was used, which permitted the recording of two gal- 
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vanometers (one attached to a thermocouple in the sample, the other in 
the standard substance) on a single record. 

In 1899, Roberts-Austen*® devised the differential thermocouple method 
for measuring differences between a material and a reference unit material. 
Burgess® elaborated the principle of the difference thermocouple and sug- 
gested simple, adaptable, and usable differential thermocouple circuits. 
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FicurE 1, Thermal reactions in clay superimposed on furnace heating rate as recorded by Le Chatelier.? 
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Ficure 2. Thermal reactions in clay 
Curve A is furnace heating curve. Curves B, C, and D show furnace heating rate with thermal reactions 
for clays superimposed. 


superimposed on furnace heating rate as recorded by Wohlin.45 


Since the work of Roberts-Austen, the difference thermocouple has been 
widely used in the field of metallurgy. Fenner,!! however, in his study of 
the stability relations of the silica minerals, reported in 1913, appears to 
have been the first to apply this differential thermocouple procedure to the 
study of thermal reactions outside of the field of metallurgy. The tech- 
nique devised by Fenner is substantially that used today, with some elab- 
oration and refinement. After Fenner’s report, the differential thermo- 
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couple method was used sporadically in the study of high-temperature 
phase changes, notably by Kracek®® and others, * and frequently in the 
study of the constituents of clay by various investigators, including Houlds- 
worth and Cobb,” Satoh,” Orcel,?4 and Kournakov and Tchernick.”8 

In 1933 and 1935, excellent papers by Orcel**: #6 emphasized clearly the 
applicability of the method to clay mineral researches. Since that time, 
many investigators in many different countries have used the method in 
clay researches. In that field, it has now become a common research tool. 

In very recent years, the method has been shown to be a valuable ana- 
lytical tool for many materials other than clays. It is still used to a con- 
siderable extent in the study of high-temperature phase changes and reac- 
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Ficure 3. A—thermocouple set-up as used today; B—thermocouple set-up as used in early work and in 
some? current work. a = Pt-10 per cent Rh wire; b = Pt wire. 


tions. It is also used widely as a plant control and raw material control 
instrument and for a quick evaluation of mineral prospect samples. 


Description of Method 


In the method as generally used today, the sample to be studied is placed 
in one hole of a specimen holder. An inert material that experiences no 
thermal reactions when heated to the temperature of the experiment, usu- 
ally calcined aluminum oxide (a@Al,Os), is placed in the other hole of the 
specimen holder. One junction of the difference couple (FIGURE 3A) is 
placed in the center of the sample and the other junction in the center of 
the inert material. The holder and thermocouple are placed in a furnace 
so controlled as to get a uniform rate of temperature increase. The tem- 
perature of the inert material increases regularly as the temperature of the 
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furnace increases. When a thermal reaction takes place in the sample, the 
temperature of the sample is greater or less than that of the inert material, 
depending on whether the reaction is exothermic or endothermic, for an 
interval of time until the reaction is over and the temperature of the sam- 
ple again comes to the temperature of the furnace. Consequently, for an 
interval of time, the temperature of one junction of the difference couple is 
different from that of the other, and an emf is set up in the differential 
thermocouple circuit which is recorded as a function of time or of the tem- 
perature of the furnace. The record may be made manually with a poten- 
tiometer or galvanometer, photographically with a galvanometer, or auto- 
matically with some electronic device. When no thermal reaction is taking 
place in the specimen, the temperature at both junctions of the difference 
couple is the same, and no emf is set up. The direction of the emf in the 
differential thermocouple circuit depends on whether the temperature of 
the sample is above or below that of inert material, and, consequently, the 
recording mechanism moves in opposite directions for exothermic and endo- 
thermic reactions. 

Ficure 4 shows a differential thermal curve for kaolinite and an equi- 
librium dehydration curve for the same material. The endothermic reac- 
tion between about 500 and 700°C. obviously corresponds to the dehydra- 
tion of the mineral. A comparison of the curves illustrates that the dif- 
ferential method is a dynamic rather than a static method. The thermal 
reactions are not instantaneous and are recorded as functions of time or as 
functions of the furnace temperature, which is continuously increasing as 
the reaction takes place. The temperature at which the dehydration begins 
corresponds to the start of the endothermic deflection. The temperature 
of the peak of the endothermic deflection varies, depending on the details 
of the procedure followed, the character of the reaction involved, and the 
material being studied. 

An attempt has been made by Speil® to analyze mathematically the dif- 
ferential method. As will be shown, certain inherent limitations, due to 
factors of experimentation and material being studied, restrict mathemati- 
cal treatment. 


Equipment Used 


Specimen Holder. A considerable variety of specimen holders has been 
used by various investigators. In early work by Le Chatelier’ and others, 
the specimen was placed in a platinum crucible. Norton,* one of the first 
to apply the method in this country, used a heavy nickel block because “it 
tends to neutralize thermal gradients.” Other investigators have used 
ceramic specimen holders, and Grimshaw'® claims they are more desirable 
than metal because the thermal reactions are larger and more sharply de- 
fined. This is due to slow transfer of heat between the test sample and 
surrounding material because of low heat conductivity of the ceramic mate- 
rial. Grimshaw!* uses a holder composed of recrystallized alumina, bonded 
with a little clay and fired to 1600°C. 

_Orcel’® and his followers® in France have placed the materials in fused 
silica crucibles. Herold and Planje?! have described a platinum cup con- 
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tainer which also serves as thermocouple, and it is claimed to simplify the 
experimental procedure, particularly for high-temperature studies. Re- 
cently, Gruver'? has claimed advantages for thin-walled platinum con- 
tainers because the high heat conductivity of the metal permits the material 
to reach furnace temperature quickly and because the thin walls provide 
little heat capacity. According to Gruver,” thickwalled containers with 
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Ficure 4, A—equilibrium dehydration curve for kaolinite; B—differential thermal curve for kaolinite 
(after Speil4?). 
higher heat capacity act as heat reservoirs which tend to lower the intensity 
and sharpness of some reactions. 

Blocks of nickel have been used by most American investigators, and, in 
general, they have been found satisfactory. Ceramics holders have been 
used widely in England. Obviously, the ceramic block must not itself react 
with the material being investigated, and for certain types of material they 
would prove unsuitable. In the writer’s laboratory, a nickel block or a 
platinum block (for high-temperature studies) has provided adequately 
sharp reactions. The metal block has the advantage of strength. This is 
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important because many materials being studied vitrify or melt in the course 
of the analysis, so that there is the problem of removing the fused button 
from the block. It would seem that a metal holder with walls as thin as is 
consistent with adequate mechanical strength might be the best type of 
block. wre s . 

Thermocouples. FicuRE 3A shows the thermocouple circuit most widely 
used at the present time. Ficure 3B shows the circuit used by Fenner" 
and Kracek2 and still used for analytical work by some European inves- 
tigators, notably the Russians.? In this latter circuit, the sample tem- 
perature rather than the furnace temperature is taken as the reference tem- 
perature. The present tendency is to use the circuit in FIGURE 3A with a 
separate thermocouple to measure furnace temperature for calibration. 
It seems probable that more reproducible results, with less variation be- 
tween results from different laboratories using slightly different equipment 
and techniques, are obtained by use of the circuit shown in FIGURE 3A. 

Early investigators used platinum-platinum 10 per cent rhodium thermo- 
couples, and such couples are very widely used today. Chromel-alumel 
couples? have been used satisfactorily for temperatures up to 1000°C. for 
many materials, and they have the great advantage over noble metal couples 
in developing higher emf and thereby providing greater sensitivity. 
Kracek?® has used gold-palladium vs. platinum-rhodium thermocouples, 
which yield about 50 microvolts per 1°C. in the elevated temperature range 
and are, therefore, about five times as sensitive as platinum vs. platinum- 
rhodium. Au-Pd vs. Pt-Rh would undoubtedly be more satisfactory for 
some materials than the sensitive base metal couples. 

Berkelheimer’ has indicated that the size of the thermocouple bead should 
be held approximately constant, because “the size of the beads in the dif- 
ference couple controls the slope of the base line of the differential tempera- 
ture curve.” 

The intensity of the thermal reactions varies greatly in different mate- 
rials. Consequently, it is desirable to have a resistance placed in series in 
the differential thermocouple circuit, which permits some variation and 
control of the sensitivity of the apparatus. 

Furnace Design. It is necessary to have a furnace with a uniform-tem- 
perature hot zone of sufficient size, so that the sample holder and any cali- 
brating thermocouple junction are all subjected to the same temperature. 
The author’ has found a horizontal furnace mounted on rollers, so that it 
can be easily rolled over an alundum tube containing the specimen holder, 
to be very satisfactory. This design permits easy access to the specimen 
holder for inserting the sample and for any thermocouple replacements. 
It also permits the use of more than one furnace with a single specimen- 
holding arm and a single recording mechanism. 


Obviously, the furnace must have enough capacity to permit the desired 


heating rate up to the desired temperature. Several authors?’ 1% 3? have 
described satisfactory furnaces in detail. 

Furnace Temperature Program Controller. A wide variety of homemade 
and automatic types of program controllers has been used to get a uni- 
form heating rate. The use of an auto transformer,’ “ driven by a motor 
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through a speed reducer in such a way as gradually and continuously to 
increase the voltage into the furnace, has been found quite satisfactory and 
economical in many laboratories for furnaces with alloy windings. The 
method does not supply a uniform heating rate for platinum-wound fur- 
naces because of the change in resistance of platinum with temperature. 
With the autotransformer-type control apparatus, heating rates varying 
less than 1°C. + at a rate of 12°C./min. can be obtained for temperatures 
up to 1000°C. 

A considerable variety’: #1: 2’. 48 of automatic controllers has been de- 
scribed. They are necessary for platinum-wound furnaces or when close 
reproducibility of heating rate is desired. In general, the automatic con- 
troller must be “tailor-made” for the furnace it is to control. This is par- 
ticularly essential for controllers of the ‘on-off type. When such con- 
trollers are used with inside-wound furnaces, or even with outside-wound 
furnaces if the furnace tube is thin, there is a tendency for the on-off im- 
pulses to get onto the differential thermocouple record, making it irregular. 

Recording Mechanism for Differential Thermocouple. In the early use of 
the method, the emf set up in the differential thermocouple was determined 
at regular intervals by means of a potentiometer, and the data were plotted 
against time or temperature. Besides being time consuming, the method 
had the disadvantage of missing some slight reactions of low sensitivity 
which might be significant. It is very important that the recording be 
continuous. 

A photographic method is used in many laboratories in which a reflecting 
galvanometer is arranged so that a spot of light is reflected from the gal- 
vanometer mirror onto a piece of photographic paper attached to a drum 
rotating at a certain speed. With this method, high sensitivity can be ob- 
tained by increasing the distances between the galvanometer and the drum. 
It has the disadvantage that the record can not be seen as it is made and 
that it requires a dark room or some rather elaborate light-tight equipment 
to handle the photographic paper. . 

Experience in the author’s laboratory on a wide variety of materials 
indicates that recording equipment must be sensitive enough to yield a 
deflection of at least two inches in length for a reaction at about 600°C. 
when the temperature difference between the junctions of the difference 
couple is no more than 10°C. For platinum-platinum rhodium thermo- 
couples, this is an emf of 0.1 M.V., and there has been difficulty getting 
recorders with this sensitivity. With base metal couples, the emf is larger, 
and automatic equipment of adequate sensitivity is available readily. Re- 
cently, Gruver,!” Kerr and Kulp,” and Rowland and Jonas” have described 
automatic recording mechanisms which have given satisfactory results. 
The automatic recorders should be equipped with variable resistances, so 
that their sensitivity can be varied and controlled. 

Calibrating Devices. Some means are necessary to determine the tem- 
perature of the thermal reactions shown by the difference couple record. 
This may be done by placing a regular thermocouple in the specimen 
block (F1cuRE 3) and taking temperature readings at specified intervals. 
In the photographic method of recording, these furnace temperatures may 
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be printed on the same paper as the difference record by placing a light in 
front of the box containing the recording drum with the light connected to 
a switch through which it can be lighted momentarily when the furnace has 
reached certain temperatures. Honeybourne” has constructed a device 
whereby the light in front of the drum is lighted by an automatic program 
controller at specified temperature intervals. 

In the photographic method of recording, it is difficult to correlate the 
height of the thermal deflections with the temperature difference at the 
time of the reactions, i.e., to set up a vertical scale for the difference curve. 
Grim and Rowland" have constructed curves applicable to their apparatus 
based on measuring the swing of the galvanometer for known temperature 
differences. 

Faust has recently presented a method of determining the temperature 
of the thermal reactions based on the inversion temperature of a and B 
quartz and the trace of the furnace heating rate recorded on the paper, 
along with the difference record, by means of a second reflecting galvanom- 
eter. When an automatic recording mechanism is used, the problem is 
that of calibrating the instrument for the thermocouples used. 


Experimental Factors Influencing Results 


It is proposed here to consider the experimental factors that can be varied 
by the operators and that can influence the results obtained. 

Packing of Sample. Gruver states that the sample needs only loose 
packing in the specimen holder. Many other users of the method, includ- 
ing the author, prefer solid packing. All samples should be packed in as 
nearly the same way as possible, so that there is the least chance of a vari- 
ation.in thermal diffusivity because of variation in density. The thermo- 
couple bead should always be of about the same size, and it is very impor- 
tant to have it centered in the center of the sample. The specimen holder 
should be placed in the furnace so that the sample and inert material are 
always in the same position. 

In the author’s experience, it is most convenient to bring the thermo- 
couples into the specimen from the side rather than from the top or bottom, 
because this permits more uniform packing, while retaining the junction in 
the center of the sample. 

Thermocouples. Relatively fine wire (gauge 22) should be used to reduce 
conductance of heat through the wire. The wire should be well insulated, 
particularly where it goes through the specimen block. 

Preparation of Sample. The particle size of the sample packed into the 
furnace is of no importance unless the material is all relatively coarse (+60 
mesh) or extremely fine-grained (—2 microns). A range of particle size 
permitting dense packing is advantageous. The author has found it good 
practice to grind all coarse samples to pass a 60-mesh sieve before analysis. 
The overall objective is to get reproducible density of packed samples in 
the block. 

Materials with the property of adsorbing water must be brought to equi- 
librium at a definite relative humidity before analysis, if reproducible re- 
sults in the initial part of the curve’are to be obtained. F1cure 5 presents 
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two curves for the same clay (prepared at different relative humidities), 
showing the difference obtained in the initial part of the curve when sam- 
ples have not been brought to equilibrium at a given RH. In studying 
clays, the author follows the practice of allowing all samples to remain in a 
desiccator with 46 per cent R.H. (saturated Ca(NO;)2-4H2O solution) for 
four days before running the analysis. 

Heating Rate. Both Norton* and Speil have investigated the effect of 
heat rate on the difference curve. As shown in FIGURE 6, the slower the 
heating rate, the broader the peak and the lower the peak temperature. 
The area under the curve for a given reaction and the temperature at the 
start of the reaction are independent of the rate of heating. 
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FIGURE 5. Differential thermal curves of same clay sample treated at different relative humidities before 
running analyses. 


The experience of many operators has shown that a heating rate in the 
range of 10° to 15°C. per minute is most satisfactory. A slower heating 
rate reduces the sharpness of the reactions and prolongs unnecessarily the 
time required for an analysis. A faster heating rate is apt to cause too 
much overlapping of reactions in samples that are mixtures. 

The heating rate must be uniform. Even small changes in the rate of 
heating, if they are relatively sudden, will be reflected in the differential 
curve. 

Furnace Atmosphere. In order to obtain reproducible results for many 
materials, the furnace atmosphere must be kept constant. For materials 
that are subject to oxidation reactions, it is important that samples be pre- 
pared and placed in the furnace in such a way that they always have the 
same ease of oxidation. Rowland and Jonas* have shown how variations 
in particle size, packing, covering of block, and dilution of sample, in addi- 
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Ficure 6. Effect of variation in heating rate on the endothermic reaction for the dehydration of kaolinite 
(after Speil42), 
tion to furnace atmosphere, influence the ease of oxidation of the iron in 
siderite and, hence, the curves obtainable from the mineral (FIGURES 7 
and 8). 
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Material Factors I. nfluencing Results 
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It is proposed here to consider the variations within the material to be 
analyzed, over which the operator has little or no control. 


CURVE "A" 


1000°C 


1000°c 


CURVE "“c" 


CURVE "D" 


100° ~—* 200° 300° 400 500 600 700° 


FicurE 7. Differential thermal curves for siderite: A—firmly packed without cover on sample holder; 
B—frmly packed without cover run in nitrogen atmosphere; C—firml: 


y packed with cover on specimen holder; 
D—very finely ground, loosely packed, and without cover (after Rowland and Jonas?°), 


800° 900°. 1000°c 


Thermal Diffusivity. The temperature at the junction of the difference 
thermocouple depends to some extent on the coefficient of thermal diffusivity 
of the material in which it is imbedded. The rate at which a point in a hot 


1042 Annals New York Academy of Sciences 

body will cool under definite surface conditions is known as thermal dif- 

fusivity” and is equivalent to an where k is the conductivity, d, the density, 
p 


and C, the specific heat. The diffusivity of the sample to be analyzed may 
be different from that of the inert material. Furthermore, the diffusivity 
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FicureE 8. Differential thermal curves for siderite show effect of dilution of sample on i 
— ue curves 0} 
containing an oxidizable component (after Rowland and Jonas‘), : oe 


of the sample may change as it is heated, because of the formation of new 
phases at higher temperatures. The effect on the differential curve will be 
to produce a sloping base line or a rather sharp offset in the curve, depend- 
ing on the suddenness of the change. In the curve for quartz (FIGURE 9), 
the base line is shifted upward after the reaction for the a and # inversion 
because of a difference in coefficient of thermal diffusivity. 
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In order to take account of this variable, Grimshaw!* would use a pre- 
calcined sample as the reference material. For materials without reversible 
heat effects this might be satisfactory. Mitchell and Henry® suggested 
making a second run on the calcined material, without disturbing the setup, 
and then subtracting the second run from that of the first. Again, any 
reversible inversions would interfere with this method. Jeffries” suggested 
using weights of sample and reference material based on their specific heats. 
This is obviously impossible when unknown mixtures are analyzed and 
would only be valuable if the sample and inert material held the same 
ratio between their specific heats throughout the temperature interval of 
the run. 


o0°c 100° 200° 300° 400° 500° 600° 700° 800° 900° 1000° 100° 


EFicure 9. Differential thermal curve for quartz. 


The matter of thermal diffusivity is very significant in the differential 
thermal method and has not yet been given the critical study that it de- 
serves. 

Changes in Volume of Sample on Heating. Shrinkage or swelling of a 

sample on heating will obviously affect the coefficient of thermal diffusivity 
and, therefore, alter the difference curve. In some materials, the sample 
may shrink away from the specimen holder with an attendant shift in the 
curve because of'a difference in heat transfer to the sample as compared to 
the inert material. 
- Variations with Particle Size. Speil® has shown that, for some materials, 
the thermal curve varies with the particle size of the component minerals 
in the size range below about 2 microns. In general (see FIGURE 10), the 
size of the thermal reaction and the temperature of the peak decrease as 
the particle size decreases. 
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In some instances, decrease in particle size is accompanied by a decrease 
in the perfection of the crystallinity of the material. Curve C of FIGURE 11 
shows a curve for a kaolin composed of well-crystallized particles of kaolinite 
of about one micron in diameter. Curves A and B of FricuRE 11 are for 
kaolins composed of particles of kaolinite less than about 0.1 micron in 
which there is some distortion of the lattice. The intensity of the reactions 
and the temperature of the peak decrease as the perfection of crystallinity 
decreases.% This influence of crystallinity may be a disturbing factor in 
some cases, but it has considerable analytical value in others. 

Variations in Quantity of Material in a Mixture. Ficure 12, after Speil,” 
presents a series of curves for varying quantities of kaolinite. The size of 
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Ficure 10 FIGURE 11 


Ficure 10. Differential therm ini i iati 
Seelle cha talter Sei). ermal curves for kaolinite showing variations resulting from differences in 


Fi Di é 
curE 11. Differential thermal curves for kaolin composed of very fine grained poorly crystalline kao- 


linite A and B, c i i 
Hate a and for kaolin composed of well-crystallized, relatively coarse-grained kaolinite C 


Wie ee and the temperature of the peak decrease as the amount of 
material decreases. In mixtures, therefore, the temperature of the peak 
cannot always be used to identify the components. 


Intimacy of Mixing of Components. Grim? has shown that the char- 


8 eee the pr i for mixtures of some minerals depend somewhat 

on the intimacy of the mixing of the components. In cl: i 

is considerable difference between th panies aan aa 

e curves for mixtures of discrete par- 

ae of clay mineral of several microns or more in diameter and an sntieaate 
erlayering of considerably smaller units, such as is frequently found in 


nature. Arbitrary curves set up on the basis of dry mixtures of clay min- 


’ ? 
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In general, the thermal reactions for the individual components decrease 
in intensity and sharpness of reaction as the intimacy of mixing increases. 

Overlapping of Reactions. Many clays contain organic material which 
burns out over a long period of time, causing an exothermic reaction which 
effectively masks any endothermic reactions which may be taking place for 
other components at the same time.!® Frequently, this difficulty can be 
overcome by firing the material beforehand to a temperature adequate to 
burn the organic material but below that of the significant reactions of the — 
other components. In some materials, the organic material can be removed 
by a solvent or H,O treatment. 

The overlapping or the concurrence at the same temperature of the 
thermal reactions characteristic of the components in a mixture may make 
the method unworkable for some materials. Grimshaw" has overcome this 


50% kaolin - 50% quartz 


25% kaolin — 75% quartz 


200 400 600 800 1,000 
TEMPERATURE, 0°C. 


Ficure 12. Differential thermal curves for kaolinite diluted with varying amounts of quartz 
(after Speil#?). 


difficulty, in some cases, by diluting the inert material with the estimated 
amount of one of the components of the mixtures, which causes the reac- 
tion for that component to be cancelled out of the difference record. Obvi- 
ously, this procedure could be used only if the identity and approximate 
abundance of the conflicting components were known. 

Influence of Minor Constituents and Impurities. The thermal curves of 
some molecules, particularly the clay minerals, are sensitive to the presence 
of very small amounts of some non-clay materials. Hendricks and his 
colleagues” have pointed out the variation in the initial endothermic reac- 
tion of montmorillonite clay minerals resulting from differences in the char- 
acter of the exchangeable base. Caillére and Hénin® and unpublished work 
in the author’s laboratory have shown that relatively large variations in the 
thermal curve for the clay minerals may result from the presence of small 
amounts of salts of various kinds which may be present in clays. 
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Quantitative Interpretations 


Care in the preparation of material for analysis and in making the analy- 
sis itself will provide reproducible curves which permit quantitative inter- 
pretations. FiGuRE 5 shows two curves, substantially alike in all details 
in the region above 250°C., which were run on the same sample at an inter- 
val of two months. In the region below 250°C. they are not alike, because 
the samples were not brought to equilibrium at the same ite 

It is obvious from the foregoing discussion of the limitations of the method 
that its accuracy for quantitative analysis is limited. The accuracy varies 
with the kind of material being analyzed, but, for many substances, values 
-+-2-5 per cent are obtainable, depending on the relative abundance of the 
components. In some mixtures, less than 5 per cent of certain components 
can be identified. 

Use of the method for quantitative analysis requires setting up standards 
by the analysis of prepared samples which contain known amounts of the 
components to be sought in the analysis. The standards must be obtained 
on the apparatus to be used in the work, and tests must be run according 
to a standard procedure. For the quantitative analysis of some natural 
materials, it is difficult to obtain satisfactory standards because of the 
difficulty in duplicating extremely fine particle size or the intimacy of mix- 
ing of the components. 

In general, the quantitative evaluations are based on the area under the 
peaks which correspond to reactions of components. It may be difficult to 
measure accurately the area of some peaks because of the difficulty of 
establishing precisely the start of a reaction or a definite base line. Meas- 
uring an area is simple for reactions which begin and end abruptly, such as 
for the a and 6 quartz transition, but, unfortunately, many reactions (such 
as those accompanying the dehydration of many hydrous micaceous min- 
erals) begin gradually, end rather gradually, and are of low intensity. 
Berkelheimer* has presented a method of measuring the area beneath de- 
flections on curves which is workable for reactions that are at least moder- 
ately sharp and intense. Dean® has used a trigonometric function of the 
peak angle for quantitative analysis which is workable also for materials 
giving relatively sharp and intense reactions. 

Recently, the International Committee for the Study of Clays (Comité 
International pour |’Etude des Argiles, CIPEA) has organized a subcom- 
mittee which will study in detail the variations in the method due to experi- 
mentation procedure, equipment, and material being analyzed. On the 
basis of this study, the subcommittee hopes to recommend a standard pro- 
cedure that will enhance the quantitative value of the method. 


A pplication of Differential Thermal Analysis 


A few differential thermal curves are presented to illustrate some of the 
various substances that have been studied by the method. No attempt is 
made to illustrate all types of the various substances that have been studied, 
nor even all the fields to which the method has been applied. 
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Ficure 13 shows curves for a mica and several clay minerals. Differen- 
tial thermal analysis has been applied more extensively to the study of clays 
than to any other material, and numerous investigators! ! 18, 3, 85, 2 have 
published many curves for specific clay minerals and natural clays. All of 
the clay minerals are characterized by an endothermic reaction somewhere 
between 500° and 700°C. due to loss of (OH) water. Some are also char- 
acterized by an initial endothermic reaction due to loss of adsorbed water, 
by a high-temperature endothermic reaction due to destruction of the lat- 
tice, and by a high-temperature exothermic reaction resulting from the 
crystallization of a new phase. The difference in the initial endothermic 
reactions of the montmorillonite clays (D and E of ricureE 13) is a conse- 
quence of a difference in the exchangeable base carried by these minerals. 
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FIGURE 13 FIGURE 14 
Ficure 13. Differential thermal curves for: (a) muscovite, coarse grained, University of Illinois collec- 
tion; (b) illite, Fithian, Illinois; (c) glauconite, Gloucester County, N, J.; (d) montmorillonite, Clay Spur, 
Wyoming; and (e) montmorillonite, Pontotoc, Mississippi. 


Ficure 14. Differential thermal curves for: (a) heulandite, Sonora, Mexico; (b) sepiolite, Ampandran- 
dava, Madagascar; (c) talc, Vermont; (d) clinochlorite, Chester, Massachusetts; and (e) wall plaster. 


For some samples, data concerning the character of the exchangeable base 
can be determined from the thermal curve, because the amount of adsorbed 
water is related to the adsorbed ion.” 

Ficure 14 presents curves for a zeolite, sepiolite, talc, and chlorite to 
illustrate that these minerals are amenable to the method. It should be 
emphasized that other types of these minerals, for example, zeolites other 
than heulandite, will give different curves. Curve E in ricurE 14 is from 
a sample of plaster and shows an initial reaction characteristic of gypsum 
and endothermic reactions at about 400°C., 500°C., and 800-900°C., result- 
ing respectively from magnesium hydrate, calcium hydrate, and a carbonate. 
The method works particularly well for hydrates and carbonates because 
they give sharp, intense, thermal reactions, and it is a simple matter to get 
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a horizontal base line. Quantitative application of the method for such 
materials is possible with relatively high accuracy. Cuthbert and Row- 
land? and Kerr and Kulp” have published curves for a variety of carbonate 
minerals. , ; 
Ficurr 15 shows curves for dolomite, limonite (ferric hydrate), bauxite, 
and alunite. The bauxite is a mixture of gibbsite which provides the endo- 
thermic peak at 300-360°C. and the halloysite-type clay mineral. The 
method works well quantitatively for bauxite (Hendricks, Goldich, and 
Nelson? have described portable differential thermal equipment used in 
bauxite prospecting). The method has the great advantage of speed, per- 
mitting a relatively close estimation in the field of the Al;Og content of 
bauxite samples in a little more than an hour. Thus, the promising samples 
can be selected quickly from the less promising or worthless ones for chemi- 


| ) 
IN 
N | \ 
| 1 | 
| 
=. 1 ‘ 
‘ue 4 t 
| ee eae hie 
it i fot iS ae 
iI ' eed Oe ee 
tl 


che or 200 300" 400 500" Gar Tor BOC BOC 1000 Noor E««OOY «OC BOC 400" SOD" 00" TOO 800r 900" 100" Noor 


Ficure 15 FIGURE 16 


Ficure 15. Differential thermal curves for: (a) dolomite, Candan, Ohio; (b) limonite, University of 
Illinois collection; (c) bauxite, Irvington, Ga.; and (d) alunite, locality unknown. 

Ficurg 16. Differential thermal curves for: (a) vitrain, diluted with 95 per cent inert (@Al:Os); (b) 
clarain, diluted with 95 per cent inert (@AlzOs); (c) pyrite, diluted with 95 per cent inert (@Al:Os); and 
(d) marcasite, diluted with 95 per cent inert (@Al2Os). 


cal analysis. Kerr and Kulp” have described equipment in which six analy- 
ses can be run simultaneously. 

Curve D in Ficure 15 is for alunite (K.0-3Al,03-4SO,-6H2O). The 
method should be applicable generally to hydrous salts of all kinds, but, so 
far, the differential thermal characteristics of such materials have been 
little investigated. 

FicureE 16 presents curves for two sulphides, pyrite and marcasite, show- 
ing distinct exothermic reactions for the oxidation of the sulphur. The 
method should be applicable generally to sulphides other than the iron 
sulphides. 

The differential thermal method is applicable also to the study of organic 
compounds. Hollings and Cobb” applied the method to the study of the 
combustion of coal. Curves A and B in FicuRE 16 show the records for 
vitrain and clarain constituents of Illinois coal made in the author’s labora- 
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tory. The exothermic reaction for such material is so great that the sample 
must be diluted with an inert material for analysis. Jordan”® has used the 
method in the study of organic materials adsorbed by the montmorillonite 
clay mineral, as shown in FIGURE 17. The amines studied by Jordan show 
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“ ae 17. Differential thermal curves of montmorillonite with adsorbed amine compounds (after Jor- 
an?6) , 


an endothermic reaction due to the fractionation of the organic molecule 
superimposed on the exothermic reaction due to oxidation. Vold* has ap- 
plied the method to the study of soaps (FIGURE 18) and has followed the 
difference-temperature record on cooling as well as heating. 

The method is applicable to the study of phase changes" taking place as 
substances are heated and to high-temperature reactions in phase equi- 
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Ficure 18. Differential thermal analysis of anhydrous crystalline calcium stearate. Curve I—base line 
obtained with Nujol in both sample and reference cell. Curve Il—heating Ca(Str)2 at 1.5°/min. Curve 
I1I—cooling Ca(Str)2 at from 2.5°/min. at 200° to 0.5°/min. at 100°. Temperatures of increase or decrease 
of AT as marked, indicate transitions (after Vold‘’). 
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librium studies, provided the change takes place abruptly and a modest 
amount of energy is involved. Ficurer 9 shows an endothermic reaction 
corresponding to the change from @ to 6 quartz, and FIGURE 11 shows an 
exothermic reaction corresponding to the formation of a new high-tempera- 
ture phase from kaolinite at about 975°C. 

Ficure 19 shows differential thermal reactions for a clay-lime mixture 
(A) and a clay-soda-lime mixture (B) made in the course of a study of a 
part of the problem of the extraction of alumina from clays. Curve A 
(kaolinite clay and calcite) shows an endothermic reaction at about 600°C. 
due to loss of (OH) water from the kaolinite and an endothermic reaction 
at 800-970°C. corresponding to loss of CO, from the lime carbonate, fol- 
lowed immediately by an exothermic reaction, which X-ray diffraction data 
indicate as resulting from the formation of gehlenite (CazAlSiO;). X-ray 
data show that pentacalcium trialuminate and dicalcium silicate have 
formed by 1400°C., but they are not reflected in the curve, which suggests 
that they formed slowly. Curve B (illite clay, soda ash, and calcite) shows 
an initial endothermic reaction due to loss of adsorbed water, a second 
endothermic reaction due to loss of (OH) from the illite, and a third endo- 
thermic reaction corresponding to loss of CO, from the carbonates, followed 
by a series of endothermic and exothermic reactions which could not be 
fully interpreted, except that the final, most intense one represents, at least 
in part, the formation of dicalcium silicate. The method has been used 
very little to study reactions of mixtures at high temperatures, but it ap- 
pears to have considerable promise, particularly when used in conjunction 
with X-ray diffraction and optical analyses to identify the phases corre- 
sponding to the thermal reactions. 
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OPTICAL PROPERTIES OF SURFACE FILMS 


By Alexandre Rothen 
Rockefeller Institute for Medical Research, New Vork, N.Y. 


“Having proved that the superficial viscosity of water was due to a 
greasy contamination, whose thickness might be much less than one mil- 
lionth of a millimeter, I too hastily concluded that films of such extra- 
ordinary tenuity were unlikely to be of optical importance.” These words 
of Rayleigh were written about sixty years ago.’ The feeling which he had 
when he first considered the optical influence of films of such tenuity is still 
the feeling of most of us when confronted with the same problem. We are 
prone to think that the yardstick should not be of a larger magnitude 
than the dimension to be measured. Influenced in part by the theory of 
resolution of optical instruments, we think of the wavelength of light as 
the order of magnitude below which no dimension can be measured by 
optical means. As we shall see, the problem of optical resolution in two 
dimensions is entirely different from that of estimating an optical path, 
which can be determined to within an extremely small fraction of one wave- 
length of light. 

In this review, we shall be concerned with optical phenomena occurring 
on reflection of light from a metallic surface coated with a film, and espe- 
cially with the application of these optical events to the determination of 
the optical thickness of the film. The optical parameters which characterize 
a film are: (1) the optical thickness, that is, the product of its real thickness 
by the index or indices of refraction; and (2) the absorption coefficient. 
Absorption occurs in metallic, that is, conductive films as well as in certain 
dielectric films whose molecules possess, for example, a system of conjugated 
double bonds. 

Let us consider the case of a transparent film deposited on a reflecting 
metallic surface. Obviously, if the film thickness is such that a path dif- 
ference of the order of \/2 occurs between the light reflected from the 
upper surface of the film and the light reflected from the metal surface, 
interference phenomena occur. On this basic fact, Blodgett and Langmuir? 
developed a method some years ago to determine the thickness of trans- 
ferred films on metallic surfaces. In their method, enough layers of a 
reference film are deposited on a polished metal slide, barium stearate, for 
example, so that the added thickness produces an optical path difference of 
slightly less than \/2 when the light is reflected at an angle of incidence of 
70° to 80°. Additional layers are then deposited stepwise so that the dif- 
ference in thickness between two adjacent steps is equal to two monolayers. 
Measurement consists in determining the angle of incidence at which two 
adjacent steps appear equally illuminated. When a film of unknown thick- 
ness is deposited on a slide covered with such an optical gauge, the angle 
of incidence must be shifted in order to reéstablish equal darkening of two 
adjacent steps. The thickness is then calculated from this shift of angle 
of incidence. The light used is linearly polarized and is vibrating perpen- 
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dicular to the plane of incidence (the reason for this choice will be dis- 
cussed later on). This method is thus based on a phase shift close to 180° 
between the electric vector corresponding to the waves reflected at the air- 
film and film-metal interphases. To be sensitive, the method requires that 
the change in phase shift versus thickness should be large, in the neighbor- 
hood of a 180° phase shift. 

This technique can be considered as a particular case of a more general 
method. Instead of producing nearly complete interference with light 


I 


II 
Fe) 


Ficure 1. The coordinates are the directions perpendicular (OS) and parallel (OP) to the plane of in- 
cidence. Ila and IIb are the ellipses obtained after metallic reflection at an angle of incidence of approxi- 
mately 70°. The ellipses are symmetrically placed with respect to the plane of incidence and 
rotate in opposite direction. 


polarized in, or perpendicular to, the plane of incidence, one can illuminate 
the surface with light polarized at 45° to the plane of incidence, that is, 
light having components both in, amd perpendicular to, the plane of in- 
cidence. In so doing, partial interference can be produced between these 
two components vibrating at right. angles, the term “interference” being 
understood in the broadest sense. 

Suppose that light linearly polarized at 45° to the plane of incidence is 
sent on a metallic surface. Before reflexion, the electric vectors can be 
looked upon as made up of two components which, in this case, are of 
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equal magnitude vibrating in phase, one in the plane of incidence (p), the 
other perpendicular to it (s). After reflection from the metal, the light 
becomes elliptically polarized. The two components suffer a phase differ- 
ence as well as a change in relative magnitude. The phenomenon is illus- 
trated in FIGURE 1, in which the plane of the figure is perpendicular to the 
direction of the light. OP and OS are the components in, and perpendicular 
to, the plane of incidence, respectively. Ficure 1, II, illustrates the condi- 
tions after reflection. OP is smaller than OS and is in advance of more 
than 90° over OS. Depending upon the quadrants in which the linear 
vibration is localized before reflection, the light becomes dextro- or levo- 
elliptically polarized. 

Let us now consider a metal surface covered with a transparent film. 
How will the film affect the ellipticity of the reflected light? It was to 
the great merit of Drude? to show that a film produces a change in phase 
and amplitude of both components ~ and s. Moreover, the two compo- 
nents are unequally affected by the presence of the film. The ratio, y, of 
their amplitude, and the difference of their phase, A = 6? — 6°, are func- 
tions of the thickness of the film. It is also important to realize that this 
property does not depend upon the anisotropy of the film. The effect 
occurs even when the indices of refraction along OP and OS are the same. 
The separate determination of a change in phase and amplitude for each 
component is not easy, except under special circumstances, as in Langmuir 
and Blodgett’s method. It is, however, very easy to determine a slight 
change in the ratio of the amplitudes of the components as well as in the 
difference of their phase, A = 6? — 6°. Qualitatively, a thin film will alter 
the ellipticity of the reflected light, but the problem solved by Drude was 
that of calculating this change in ellipticity as a function of the thickness 
and index of refraction of the film. 

Starting from Maxwell equations and the appropriate boundary condi- 
tions, Drude obtained the following equation: 

D,8 ps p,8 —iz 
ete rar S: (1) 
eh ees ae 


In this equation, I and p are the amplitudes of the electric vectors (square 
root of light intensity) of incident (J) and reflected (R) light; 6 is the phase, 
rz is the real Fresnel coefficient for reflection at the film surface, and fm is 
the complex Fresnel coefficient for reflection at the metal surface. The 
superscripts, ?*, indicate that this formula holds for each component sepa- _ 


rately. The variable x is equal to 4rn cos eS, in which 7 is the index of 


refraction of the film, y is the angle at which the light is reflected at the 


film metal interphase, / is the thickness of the film, and ) is the wave length 
of the light used. 


Writing the above formula for the two components, # and s, separately, 
and taking the ratio, one obtains: 
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In this equation, y = ne and A = 6? — 6°, 
p? 


It is apparent from the above formula that y and A are periodic functions 
of the thickness / of the film. It is difficult to say more than this qualita- 
tive statement without expanding the equation. The equation was devel- 
oped by Drude using the first term of a serial development. He obtained 
the following linear relationships between A or y and the film thickness: 
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In this equation, A’ and y’ are the values of A and y, when there is no 
film; A, B are functions of the angle of incidence and the optical constants 
of the metal; and C = x? cos? g, n being the refractive index of the film 
and ¢ the angle of incidence. Tronstad‘ was the first to test such relation- 
ships with thin films of known thickness. He used films of fatty acids 
adsorbed on a mercury surface and found good agreement between experi- 
mental and theoretical values. With the advent of Langmuir and Blodgett’s 
method for the transfer of films by the dipping method, it became much 
easier to test Drude’s formulas. We used films of multilayers of barium 
stearate deposited on stainless steel slides. The following values could be 
calculated from formulas (3) for one monolayer of barium stearate 24.4 A 
thick and in which 75g93 = 1.495 (TABLE 1). 


TABLE 1 
Angle of incidence: 49.2 69.2 74.2 
y—-y= 0.00 0.34 0.61 
A’-A= Ley 3.6 3255 


It was found® that these values are in complete agreement with the ex- 
perimental data obtained for thin films. At large angles of incidence, there 
is an appreciable deviation for films thicker than 50 A, although the agree- 
ment is still good at 50° incidence with films up to 250 A in thickness. 

Drude’s approximations, or similar approximations,® are thus valid for 
films which are thin with respect to the wavelength of light, but the ap- 
proximations obviously cannot hold for thick films. Hauschild’ obtained a 
better approximation by using a two-term development. This procedure, 
however, is unsatisfactory. Fortunately, formulas (1) or (2) can be developed 
explicitly without approximation. This has been done by Vasicek* for the 
case of films deposited on a transparent medium. The same method can 
be used for films deposited on metallic surfaces,® but the algebraic equations 
are a little more cumbersome since the coefficients r”,, are complex and 
must be split into two parts 77%, = X?*-++ iY2*. The development of 
formula (1) consists of separating the real from the imaginary quantities. 
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Writing this formula for the two components separately, and substituting the 
subscripts ; and 2 for the superscript ”, and subscripts 3 and 4 for the super- 
script %, one obtains: 
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Taking the ratio, 
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In this last equation, 
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B= -xX*{1 — r"] sina + Y?[1 — r| COS 2; 
and B’ = X*[1 — rs] sna — Y*[i — rs | Cos &. 


The coefficients r?> can be calculated from Fresnel equations, since the 
angle of incidence and the index of refraction of the film are known. The 
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quantities X?*, Y?.* are related to the phase shifts 5m, 5m (abbreviated 5”’*) 
of the components after reflection from the uncoated metal surface by the 
following formulas: ig’?* = Yr#/XpP.s, ph? — [Xpe 4 ye}, The correct 
sign should be taken for the coefficients X?*, Y?*, It turns out that X? 
is positive, Y? is negative, X* is negative, and Y® is positive. The quan- 
tities 67;° and A, = 6?, — 6%, as well as p?'* can be computed from the optical 
constants of the clean metal surface by means of equations which can be 
found, for instance, in the Handbuch der Physik, 20: 242. Thus, fg A can 
be calculated for any film thickness from equation (4), since the coefficients 
A, A’, B, and B’ can be computed. Ficure 2 shows that there is good 
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Ficure 2. Phase shift A produced by films of Ba stearate at 50°, 69.2°, and 80° incidence. The dotted 
curves are theoretical, and the solid and dashed lines are experimental. 


agreement between the experimental and theoretical curves of A plotted 
against the thickness of film. As indicated in (4), the change in phase A is 
not a sinusoidal function of x, since the function contains sin?x, cos’, sin x 
cos x, terms corresponding to a periodicity of 2x. The theoretical number 
of monolayers of barium stearate which produces a path difference, of 5893 
A at the angles of incidence of 50°, 69.2°, and 80° are 93.3, 102.5, and 105, 
respectively. These figures compare favorably with the experimental values 
of 92.5, 102.2, and 105 monolayers at which the periods describing the 
change of phase are completed. Furthermore, the curves of FIGURE 2 show 
that the change of phase with increasing thickness is especially large for 
films equivalent to 40 to 50 monolayers of barium stearate. In that range 
of thickness, an increment of 24 A brings about a change of phase A = 48° 
at 80° incidence. Since a change of phase A can be measured within +0.04°, 
it follows that a variation in thickness of less than 0.1 A can still be appre- 


1060 Annals New York Academy of Sciences 


ciated. When the film is about 600 A or 1500 A thick, the phase changes 
very little with the small increment in thickness. 
It is interesting to note that, because of the skewness of the curve, the 


coefficient a which is maximum at 45 layers, is also large at 95 layers, 
6 


but appreciably smaller than at 45 layers. This means that 6? and 6° taken 
separately are different functions of x. These absolute phase shifts 6? and 
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Ficure 3. Absolute phase shift 5”, 5° of each component produced by films of Ba stearate at 69.2° incidence. 


6° can be calculated from equation 1 as A was calculated above. The cal- 
culations give: 


a102 + Bi Be aaa + B3Bs- 
The results have been plotted in ricure 3 for the angle of incidence of 
69.2°. The difference of the two curves gives the corresponding A curve 
seen in FIGURE 2. The curves of FIGURE 3 show that there is a phase shift 
of 180° for the components s and 4, for thicknesses equivalent to 46.8 and 
94.1 monolayers of barium stearate, respectively. This means that a mini- 
mum does occur in the intensity of the light corresponding to the compo- 
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nents s and p at these thicknesses. The figure also shows that the curve 8° 
is much steeper than the curve 6” when their respective phase shift is 180°. 
This is why it is much better to use the component s, rather than p, in the 
method of Blodgett and Langmuir described above. The required phase 
change of 180° occurs for both of them at the appropriate thicknesses, but 
the sensitivity (change in 6 versus change in thickness) is much greater for 


Ficure 4. Ellipticity of reflected light at 50° incidence produced by films of Ba stearate from 0A up to 
the equivalent of 17 wavelengths of light (45893). Arrows indicate beginning of the first cycle. All cycles 
are described clockwise. The points indicated on the initial cycle correspond to thickness intervals of 146 
A (six monolayers of Ba stearate). 


s. Blodgett and Langmuir? indeed have noticed that the useful range of 
angles of incidence is much smaller for p thans in accord with these theo- 
retical curves. 

As already indicated, there is a corresponding periodic change in the 
value of y. Since the change in y is, in general, smaller than the change 
in A, we shall not consider it here. 

Films Producing an Optical Path Difference Greater Than One Wavelength 
of Light. With the increase in the thickness of the film, formula (2) predicts 
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that the changes in y and A are truly periodical and repeat the same cycle 
each time the thickness is equivalent to a path difference of one wavelength 
of light. To test this point, experimental data were obtained for film 
thickness up to 40,000 A of barium stearate. Both w and A were deter- 
mined. The results obtained at 50° are represented graphically in FIGURE 4 
according to the method of Poincaré. The plane of the figure is the com- 
plex plane. The coordinates u + iv of each single point are: 


u = cot ycos A; v= cotysin A; and u + iv = cot pe . 


With this representation, all points on the axis correspond to linear vi- 
brations, whereas all points on the v axis correspond to ellipses whose axes 
are oriented along OS and OP. The points, u = 0, and v = +1, are two 
circular vibrations of opposite sense. It is apparent that each cycle corre- 
sponds to an optical thickness of one wavelength of light, but the successive 
cycles precess clockwise around the origin as the approximate center of rota- 
tion. This lack of superposition cannot be predicted by formula (2). It 
may be due in part to the birefringence of the film which has hitherto been 
neglected. 

We have reviewed so far the optical phenomena occurring with trans- 
parent films. In the case of absorbing films, the general formula (1) of Drude 
is still valid except that the exponential coefficient x includes the absorption 
coefficient. This case will not be discussed since no experimental data have 
yet been obtained for absorbing films of known thickness. 

Apparatus to Measure the Ellipticity of Reflected Light. The basic prin- 
ciple of all the apparatus developed up to the present depends on a change 
from the elliptically polarized light to linearly polarized light by the intro- 
duction of a supplementary phase shift by means of a compensator (Babinet 
compensator or mica \/4 plate). In order to be sensitive, the apparatus 
should be equipped with a half shadow device. ._ In our own apparatus, built 
specifically for thickness determination, the half shadow is produced on the 
reflecting surface itself.° With such an apparatus it is an easy matter to 
determine an optical thickness of a thin film to lessthan1 A. The theoretical 


limit is -£0.1 A under the most favorable conditions (when = is max- 
imum). 

Conclusions. It has been shown that a film as thin as 0.1 Angstrom unit 
produces some noticeable effect on the parameters characterizing the light 
reflected from the surface upon which the film is deposited. The thickness 
of the film can be calculated from the magnitude of this optical effect. It 
may be well to emphasize that such a thickness is an overall thickness. 
Apparatus built with the purpose of determining such small thickness may 
be compared to a one-dimensional microscope, with a one-dimensional 
resolution far greater than that achieved in the best electronmicroscope. 
Unfortunately, it is only one dimensional. The resolution in the other 


two dimensions perpendicular to the thickness is limited by the ordinary 
laws of optics. 
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* In this article, a misprint occurred in the Drude formula appearing on page 67. The formula should 
read A’ — A, instead of A — A’. 


MEASUREMENT OF PROPERTIES OF THIN FILMS ON 
CHROMIUM BY THE REFLECTION OF 
POLARIZED LIGHT 


By J. B. Bateman and Margaret W. Harris 
Camp Detrick, Frederick, Maryland 


Thin films coating a plane polished metal surface can be characterized 
in terms of their effect upon the optical properties of the surface. Drude’ 
demonstrated sixty years ago that the effect is readily measurable and the 
experimental arrangement used by him in his studies of metallic reflection 
has survived without substantial modification through the several subse- 
quent investigations in which optical examination of surface films has been 
attempted. If a specimen of the clean polished metal is available, the meas- 
urement generally presents no great difficulty, and it may afford a con- 
venient method of following surface changes brought about by various 
agencies. More often, it is desirable to derive values for the thickness and 
refractive index of the film from the measured quantities and to interpret 
these ultimately in terms of molecular dimensions and molecular orientation. 

Difficulties of interpretation have been the chief drawback in the wider 
use of polarization optical methods in the study of surface configurations. 
Several approximations of the general equations, applicable only to very 
thin films, are available. These should describe the observed behavior 
of such films, provided it is possible to specify the boundary conditions 
which operate. Until very recently, however, no experimental test had 
been made, because acceptable independent methods of determining the film 
constants were not available. 

Increasing realization of the importance of thin films in electrochemistry, 
electronics, applied optics, and biophysics has made it desirable to establish 
the range of application of the theoretical equations’ of the optical study 
of films. Fortunately, the circumstances which give particular interest to 
the problem have also provided an excellent experimental method with which 
to attack it. Blodgett? has shown how to transfer almost any desired num- 
ber of monolayers of stearic acid and other long chain fatty acids from the 
liquid surface on which they were formed to a glass or metal plate, and has 
measured the thickness and refractive index of the transferred films.2:5 A 
study of their effects upon the optical properties of the metal base, with 
an obvious bearing upon the applicability of the theoretical equations, 
would also clearly serve as an empirical means of correlating the observed 
effects with film thickness and refractive index in the event that discrepan- 
cies between experimental and theoretical data should render the equations 
of doubtful practical value in their present form. 

This paper presents the results of polarization measurements on aliphatic 
acid layers transferred to chromium by Blodgett’s method. The data were 
obtained with the above considerations in mind, and with the special object 
of determining the conditions most favorable to the precise simultaneous 
measurement of thickness and refractive index. To this end, it proved 
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necessary to extend the measurements to thicker films than those which 
had been used for the purpose of comparison with the approximate theo- 
retical equations, most of which apply only when the film thickness is much 
less than the wavelength of the light in which they are observed. 


Approximate Equations of Film Optics 


Polarimetric methods permit the measurement of the relative changes in 
phase and amplitude brought about by reflection of a plane-polarized beam 
resolved into components whose electric vectors vibrate parallel (ray op) 
and perpendicular (0s) respectively to the plane of incidence. If the two 
components are initially equal in phase and amplitude, the ellipse defining 
the reflected vibrations is represented by the relative reflection coefficient, 
tan y-exp (7A), the real part of which is the ratio of the reflected ampli- 
tudes, R,/R;,, while the imaginary part, A, gives the difference between 
the absolute phase changes, (5, — 4;). 

The effect of a thin film upon the relative reflection coefficient of a metal 
is conveniently represented by 


tan p 
tan Wo 


exp i(A — Ao) 


in which the subscript 0 refers to the uncoated metal. Typically, when 
light is reflected at an angle of incidence 70° from the surface of a metal 
covered with a film of refractive index 1.5 and thickness 20 AU, tan W/tan Wo 
is about 1.01 and (A — Ao) has a negative value of two or three degrees. 
The available theoretical equations predict correctly the sign and order of 
-magnitude of these changes but differ among themselves as to the exact 
relationship with film thickness and refractive index. 

The experimental studies reported here have been guided by reference to 
only two of the theoretical approximations, the earliest by Drude® and a 
recent modification given by Lucy.'® 

Drude’s equations are linear in film thickness, L: 


2 
BEE seed PeMP Pate 5 : (1 — 2’ cos’ ¢) (1) 
tan Yo nN 
n —1 

A — Ay = —BL a (cos’p:-— a) . (2) 
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In these equations, g represents angle of incidence, 2 is the wavelength, 
m is the film refractive index, and mp and xy are the refractive index and 
absorption coefficient of the metal base, respectively. The presence of the 
term 7? cos? g in equation 1 establishes the possibility, in principle, of deter- 
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mining both Z and from a single measurement of tan y//tan yo and A — Ao, 
these quantities being related by the equation: 


cos? g — ay tan y me | 
tan Wo 


—(A — Ay) = (3) 


ai(1 — n’ cos’ ¢) 


FILMS ON CHROMIUM 
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It is apparent from equation 3 that (1 — n? cos * g) is least sensitive to 
changes of film refractive index when the latter approaches unity; the lines 
in FIGURE 1 illustrate this and show that the method is rather unsuitable 


for measurement of the refractive index of films of unknown thickness in 
the range below 1.5. 
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Equations 17 to 22 given in Lucy’s paper!’ can be rearranged in the fol- 
lowing form for the convenient calculation of theoretical curves: 


tan” Wisoe 2 2 
raters api burt (4) 
cos (A — Ao) = (1 — My) [1 — M)? + M2}? (5) 


In the equations, 


4r singtang n —1 


M, = — -—*—* . —__. -L .[(as — Br)n’ + (6° — a’)s + 2a6r); 
% 4r singtang nm = 
N r2 + 32 ne 


a = nl — x); B = 2nimo;r = (a — 1)(a — tan? y) — 6; and 
s = B(2a — tan? » — 1). 
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As in the Drude equations, the film thickness term always involves the 
product L-(n? — 1)/n?, while the refractive index also occurs elsewhere in 
the definitions of M; and M;. It is also evident that the equations can 
be solved in principle for ~ and ZL, and that the refractive index effect is 
not large in the case of films of refractive index less than 1.5 (see FIGURE 1). 
The curves of FIGURE 1 also show that, in the case of films of low refractive 
index, although the Drude and the Lucy equations give concordant values 
of the limiting slope at zero film thickness, deviations from the linear relation- 
ship may be expected to appear in the case of the phase change when the 
film thickness exceeds 50 AU, and in the case of the amplitude factor at 
thicknesses as low as 10 AU. When the film refractive index is greater 
than 2.5, discrepancies are already evident in the values for the limiting 
slopes. 


Experimental Technique 


In outline, the following experimental procedure has been followed: 
Preliminary observations establish the optical constants of a clean metal 
surface. The desired number of aliphatic acid monolayers is transferred 
to the surface by dipping. The polarimetric observations are repeated. 
Further monolayers are transferred and the process continued until the 
desired film thickness is attained. Detailed descriptions of each opera- 
tion are available in the literature and the following notes refer to these 
accounts and to the modifications that have been made in the course of the 
present work. 

Metal Slides. In the first experiments, optically flat polished stainless 
steel plates formed the reflecting base. It was soon found that the plates 
could not be polished in such a manner as to give uniform values of the 
optical constants. Furthermore, the constants were altered by the process 
of conditioning the plates for the uptake of an initial film when raised 
through an aqueous interface occupied by a monolayer of barium stearate. 
As a result, the apparent effect of the first one or two stearate monolayers 
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was highly erratic, since the optical constants of the conditioned surface 
could not be measured. These difficulties were not encountered when glass 
slides, coated with chromium by condensation in vacuo, were used instead 
of steel. Satisfactory results were obtained with ordinary lantern slides 
cut to the desired size. After cleaning in bichromate, followed by rinsing 
in distilled water, the slides were heated in air or in vacuo for several minutes 
and were then exposed to chromium vapor. A strongly adherent layer was 
deposited. The optical properties of the surfaces were more reproducible 
than those of polished steel. Furthermore, the bare chromium will, with 
few exceptions, pick up one fatty acid monolayer when lowered through a 
liquid interface covered with such a film, and a second when raised. The 
optical effects of this Y-layer (Blodgett?) were reproducible and consistent 
with those of subsequent Y-layers, and were brought about wholly by the 
transferred monolayers, not by mere immersion of the slide in the substrate. 
Occasionally, freshly prepared chromium plated slides are wettable. In 
this condition, they pick up a monolayer when lowered through the film- 
covered interface and release it when raised. Within a few hours, however, 
they acquire hydrophobic properties without undergoing any measurable 
change in optical constants. 

Transferred Monolayers: Skeleton Films. Barium stearate-stearic acid and 
barium palmitate-palmitic acid films were made by spreading the respective 
acids dissolved in benzene upon water containing barium chloride and buf- 
fered with bicarbonate at pH 7 (solution “a” of Blodgett and Langmuir‘). 
Cadmium arachidate layers were prepared at pH 5.7 as described by Blodg- 
ett. After observations had been made upon the continuous films of cad- 
mium arachidate, they were skeletonized by soaking for two minutes in 
acetone. When measurements on the skeleton films were complete, the 
original refractive index was re-established by treatment of the slides with a 
drop of Nujol (bibliographical reference 5, page 980), and a further series of 
measurements provided assurance that no significant loss of thickness had 
accompanied the skeletonization. 

Polarimetric Measurements. The polarimetric assembly consisted of the 
usual elements mounted upon a spectrometer base in the following order: 
sodium arc, filter, iris diaphragm, collimator pinhole, collimating lens, Glan- 
Thompson polarizer, reflecting plate, Soleil-Babinet compensator, Glan- 
Thompson analyzer, focusing lens, and eyepiece pinhole. In the earlier 
measurements, Rothen’s half-shade device was employed," consisting of an 
optical gauge of barium stearate* incorporated in the reflecting base. This 
method gave very satisfactory empirical calibrations in terms of film thick- 
ness. Matching of the fields occurs, not at some unique setting of compen- 
sator and analyzer azimuths, but over a wide range of settings. Conse- 
quently, it is possible to select arbitrarily some constant setting of the 
compensator and to determine the variation of the analyzer readings at which 
the fields match, as a function of film thickness. This procedure is less 
suited to the exact determination of the phase and amplitude factors. 


— Cae excellent half shade can also be made by condensing a thin layer of silicon dioxide on one half of 
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Rothen” achieved the latter end in his studies of thick films by adjusting 
both quarter-wave plate and analyzer to “a position of maximum uniform 
darkness.” In our hands, this procedure proved insensitive. Better results 
were obtained by discarding the optical gauge and adjusting compensator 
and analyzer to extinguish the reflected beam. Careful alignment of the 
optical elements of the instrument gave scale correction data, and further 
allowance for the effects of any remaining lack of symmetry was made by 
taking readings with polarizer settings of 45° and 315° (Drude). The 
equations for calculating tan y and A from the observed compensating azi- 
muths are given in convenient form by Lucy.! 


Results and Discussion 


Optical Properties of Chromium. The optical constants of uncoated chro- 
mium were calculated from A and tan y by means of the approximate for- 
mulae (bibliographical reference 7, page 544, and bibliographical reference 9, 
page 363). Theoretical curves for comparison with our film data were ob- 
tained by using the values, m) = 3.10, and xo = 1.22, as the optical constants 
of the uncoated chromium base. When the calculations were made, these 
were selected as the most frequently occurring values. Later averages, 
based upon measurements of 48 plates, were 7 = 2.963 and x = 1.132. 
On the other hand, calculations of A and tan y for filmed plates were made 
using the appropriate values of mo and xo for each individual plate. The 
extreme values, approached only in a few instances, were m = 2.626, 3.264; 
and xo = 0.986, 1.264. 

Behavior of Thin Aliphatic Acid Films on Chromium. An example of 
the observed optical behaviour of barium stearate films on a single chro- 
mium-plated slide, as the film thickness is increased in steps of two mono- 
layers from zero to one thousand AU, is shown in FIGURE 2. The observa- 
tions were made in sodium light at an angle of incidence of 70°. Whereas 
A reaches a maximum value around 600 AU and tan y begins to increase 
sharply at about 800 AU, the approximate theoretical equations predict a 
continuous decrease of A and a steady increase of tan w (FIGURE 1). The 
comparison of experimental and theoretical data must be confined to films 
less than 400 AU thick (optical thickness of the order 0.1 2). 

In FicurE 3, the experimental data for barium stearate and barium pal- 
mitate films are plotted, together with the appropriate theoretical lines calcu- 
lated from the Drude and Lucy equations. The scale of thickness, L, for 
the experimental data is based upon values of 48.8 and 44.0 AU, respectively, 
for the thickness of double molecular layers® of barium stearate and barium 
palmitate, while the calculated curves are for films of a refractive index 
of 1.5. 

Both the Drude and Lucy formulations are consistent with the experi- 
mental values of A — Ay for film thicknesses up to 100 AU. At this point, 
deviations from the linear relationship become evident; up to 200 AU, the 
points follow the Lucy curve within the limits of experimental error. De- 
viations from linearity in the case of tan y/tan yo are anticipated by the 
Lucy formulation for quite thin films, and, in practice, are already suggested 
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by the points for one double layer of barium stearate of a thickness of 
48.8 AU, and are quite definite at 88 AU. Quantitatively, it seems that 
the Lucy formulation for tan y/tan wo over-corrects for deviations from 
the linear relationship. 
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tematic differences between the slides and partly of irregularities in the 
effects of successive double layers. The former are partly accounted for 
by differences in the optical constants of the uncoated metal surfaces. 
Since such differences are of importance, the data should be compared with 
theoretical zones corresponding to the range of values of mo and xo encoun- 
tered in the measurements, rather than with the single curves shown in 
FIGURE 3. These differences are less important, however, than the variable 
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Ficurst 3. Thin films of barium stearate and barium palmitate on chromium: comparison of experi- 
mental and theoretical data for wavelength 5893 AU, angle of incidence 70°. Lines calculated for films of 
refractive index 1.5 upon a reflecting base with optical constants m, 3.10; xo, 1.22 from equations of Drude 
and Lucy. Experimental points are for films of barium stearate and barium palmitate, each type of symbol 
representing a single reflecting plate. The optical constants of the metal base had the following values for 
the plates used in these measurements: barium stearate plates, m = 3.130, 3.137, 3.177, 3.264; xu = 1.167, 
1.168, 1.154, 1.139; barium palmitate plates, 7 = 3.119, 3.080; xo = 1.164, 1.208. 


effects of successive double layers upon a single slide. Irregularities in the 
character of the monolayers when spread on the original liquid substrate, 
and variations in the configuration they assume when transferred to the 
chromium slide, are probably responsible for the most serious optical incon- 
sistencies, the causes of which have not been fully determined. Variations 
in the purity of the water used in the aqueous substrate, in the rate and 
temperature of transfer, and in the dust content of the air are contributory 
factors which will require careful control in any attempt to obtain more 


precise data. 
Errors referable to variations in the metal base and the transferred films 
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outweigh those arising from the technique of optical measurement and ras 
uncertainties in the values chosen for film thickness and refractive index. 
The thickness scale may require slight revision. The value 48.8 AU. for 
the thickness of a double molecular layer of barium stearate was obtained 
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Ficure 4. Effect of refractive index difference on optical behavior of thin films on chromium: com- 
parison of experimental and theoretical data for wavelength 5893 AU, angle of incidence 70°. Continuous 
lines are calculated from Drude equations for films of refractive index 1.5 and 1.23. Broken lines, from Lucy 
equations. Small experimental points are for mixed films of cadmium arachidate-arachidic acid; large points 
are for the cadmium arachidate films prepared by skeletonizing the mixed films with acetone. \ Each individ- 
ual small circle represents a separate plate (contrast data in FIGURE 3). 


by Blodgett and Langmuir‘ by observation of the intensity minima of the 
ordinary ray, R, , reflected from barium stearate films on chromium-plated 
brass, assuming the refractive index of the films for the R, ray to be 1.491, 
the value measured for similar films on glass. The length of the double 
molecule of stearic acid in monoclinic crystals is 49.6 and that of rhombic 
sodium stearate, 51.8 AU.® Although the effect of the barium atom in the 
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mixed stearic acid-barium stearate films is not known, these figures suggest 
that the chains in the transferred layers on chromium are not quite vertical. 
The tilt has been taken into account in the values 44.0 AU assumed for the 
double molecule of barium palmitate and 53.6 (q.v.) for cadmium arachidate. 
These figures are subject to slight revision when more direct measurements 
of thickness become available. The refractive index, 1.5, used in calcu- 
lating the theoretical curves, is the effective value for the angle of incidence, 
70°, calculated from m, and ng, the refractive indices of the birefringent 
film as measured by Blodgett and Langmuir.® 

Skeleton Films of Cadmium Arachidate. Films of cadmium arachidate- 
arachidic acid on chromium showed the same behavior in relation to the 
theoretical curves as barium stearate and barium palmitate (FIGURE 4), 
although the films were more erratic in their optical properties. When the 
mixed films were converted by the action of a solvent to skeleton films of 
the metallic arachidate, A — Ao and tan W/tan Wp changed to values charac- 
teristic of much thinner continuous films. This effect is predicted by the 
theoretical equations. The curves for films of refractive index 1.23, the figure 
given by Blodgett* for skeleton films prepared by her method, are shown in 
FIGURE 4 for comparison with the experimental points. The changes are 
somewhat greater than those anticipated on theoretical grounds, and are 
more marked in the case of the thinner films. 

Treatment of skeleton films with Nujol restored their optical properties 
to those of the original continuous films, except in the case of plates covered 
with only two monolayers, which seemed to be removed completely by the 
briefest contact with acetone. With this exception, therefore, it seems that 
the exaggerated effects of the solvent upon the thinner films cannot be 
explained by detachment of entire monolayers from the upper surface, but 
must be due to removal of cadmium arachidate or residual amounts of the 
free acid from the body of the film. Observations on thicker films sug- 
gest, however, that this process does not occur uniformly throughout the 
depth of the film, since prolonged immersion in solvent does result in a 
decrease in the measured thickness of the Nujol-filled skeleton (bibliographi- 
cal reference 5, page 980). 

A Semi-Empirical Linear Relationship Based on Lucy’s Equations. Lucy’s 
equations lead to a simple expression involving (A — Ao) and tan y/tan Yo 
as a linear function of Z and (n? — 1): 


[cos (A — Ac) + ’ sin (A — Ap)]-tan y/tan po = Rk’ L(n? — 1) —- 1 (6) 
In this equation, ; 
= [(? — B)s + 2aBr]/[(6? — o?)r — 2aBs] (7) 


pr = sin g tan g-B(G + 8) /[(6" — or — 2a (8) 


The symbols have the definitions given for equation 5. Under the condi- 
tions of our experiments, k’, by equation 7, should be 0.354. A plot of the 
data for barium stearate, using this value, gives the upper curve in FIGURE 5. 
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The value of k’ can be adjusted empirically to give a straight line for films 
up to 400 AU thick. The best value of k’ is unity, the slope of the straight 
line, k’, being about one order of magnitude larger than that defined by 
equation 8. 

"The same line fits the data for cadmium arachidate (FIGURE 6) and the 
same value of k’ (unity) is applicable to the skeleton films, as shown by the 
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FicurE 5. Test of equation 6 based upon Lucy’s approximation: Upper curve: values of f(y, A), using 
k, calculated from the theoretical equation. Lower curve: empirical version of equation 6 using k’ = 1.0. 
The experimental data are the same as those plotted in FIGURE 3. 


upper line in FIGURE 6, which is consistent with the experimental points up 
to 800 AU. Comparing the slope of this line with that obtained for the 
continuous films of cadmium arachidate-arachidic acid, 1.25 is derived as 
the refractive index of the skeleton films, in good agreement with the value 
of 1.23, calculated by Blodgett? for films from which one half of the substance 
has been dissolved. 

Measurement of Thickness and Refractive Index of Thin Films. For films 
of known refractive index around 1.5, measurement of (A — Ao), followed 
by application of Drude’s or Lucy’s approximate equations, will permit good 
estimates of film thickness, if less than 100 AU. When it is between 100 
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and 200 AU, Lucy’s equation is more reliable. In the particular case under 
consideration (n = 1.5), greater accuracy can be achieved empirically, 
since it is easy to carry out a calibration with barium stearate layers. It is 
then unnecessary to determine absolute extinction positions or to calculate 
A — Ao. The elegant half-shade method developed by Rothen" can be 
used. In other instances, when the refractive index of the film can be de- 
termined by independent methods and is known to differ appreciably from 
1.5, the range of usefulness of the two equations is not immediately ap- 
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Ficure 6. Test of the empirical version of equation 6 with experimental data for mixed continuous films 
of arachidic acid-cadmium arachidate (light circles) and for films of cadmium arachidate produced by skele- 
tonizing the mixed films with acetone (heavy circles). In each case, the value k’ = 1.0 was used in calcu- 
lating f(~, A). The lower straight line is the same as that shown in the lower portion of ricuRE 5. The 
upper line was drawn by inspection through the heavy circles. The experimental data are the same as those 
plotted in FIGURE 4. 


parent from the present measurements. It can be safely said, however, (see 
FIGURE 1) that both equations can be applied to increasingly thick films as 
the refractive index is decreased. 

Measurements of tan ¥//tan yp are of little use in determining film thick- 
ness by use of the theoretical equations. The failure of the equations to 
describe correctly the course of change of tan ¥/tan y implies also that 
they cannot be used for the simultaneous calculation of film thickness and 
refractive index, nor are empirical methods of much value in this connection, 
since the relationship between A — A» and tan y/tan Yo is very insensitive 
to changes of refractive index in the range from 1.0 to 1.5. This was sug- 
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gested by the curves in FIGURE 1 and is evident in the data plotted in the 
tom section of FIGURE 4. 
is some applications of the optical study of thin films, the film thickness 
is known or can be determined by an independent method. The success of 
equation 6, in its empirical form, in describing the behavior of both con- 
tinuous and skeletonized films, suggests the possibility that this equation 
will be useful in measurements of refractive index. We hope to test its 
validity for condensed films of highly refracting substances. ae 
Thick Films of Barium Stearate (up to 1000 AU). The progressive in- 
crease in ellipticity of the light reflected from filmed chromium plates with 
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Ficurre 7. Thick films of barium stearate on chromium: ¢f. bottom section of FIGURE 3, which is identi- 
cal to the initial part of this figure. Here the results of extending the observations to film thickness 900 
AU are plotted. The continuous line was drawn through the experimental points by inspection. Film 
thickness is indicated by numerals on diagram. See text. 


increasing film thickness is illustrated in FIGURE 2, together with the changes 
in the compensator and analyzer azimuth settings at which the ellipses are 
compensated, and the corresponding values of A and tan y. The curves 
draw attention to the fact that the sensitivity of the azimuth readings to 
small increments in thickness increases gradually up to a total thickness of 
about 300 AU, remains roughly constant thereafter, until, at approximately 
700 AU, it starts increasing sharply toward the point (800 AU) at which the 
reflected light becomes circularly polarized. The readings are then indeter- 
minate, and, in a zone on either side of this point, the high sensitivity is offset 
by the poor definition of the intensity minima, so that the azimuth readings 
are subject to considerable error. Beyond this zone, for films more than 
900 AU thick, is a region (see also FIGURE 7) in which great sensitivity of 
the phase and amplitude parameters to small changes in thickness is asso- 
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ciated with conditions which enable precise adjustments of the instrument 
to be made. Rothen" has indicated this region as one in which thickness 
increments of 0.1 AU might easily be detected. 

Thick Skeleton Films: Measurement of Refractive Index of Thick Films. 
The minimum value of A occurs in the case of barium stearate films at a 
thickness of about 600 AU and is about 80°, corresponding to —50° for 
A — Ao. At this point, tan y/tan yo is about 1.65. Since the rise of 
tan y/tan Yo is continuous, the relationship between the two parameters 
takes the form shown in FicuRE 7. Essentially the same result was ob- 
tained for cadmium arachidate (FIGURE 8). When the arachidate films 
were skeletonized, it was found that the minimum value of A — Ao was 
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Ficure 8. Effect of film refractive index on optical behavior of thick films on chromium. The initial 


portion of this diagram is the same as the bottom part of ricurr 4. Here the results of extending the ob- 
servations to thicker films are plotted. Film thickness is indicated by numerals on diagram. See text. 


only —40°. It occurred when the film thickness was about 800 AU and 
tan y/tan wo was 1.4. Up to this point, the curves of A — Apo against 
tan w/tan Wo for the continuous and the skeleton films had been indis- 
tinguishable within the limits of experimental error; thereafter, they di- 
verged (FIGURE 8). Thus, the regions of great sensitivity to changes in 
film thickness are also characterized by a greatly increased dependence of 
the phase-amplitude relationship upon film refractive index. A single ex- 
perimental observation upon a film represented by a point situated any- 
where in the upper third of rrcuRE 8 should therefore be sufficient to define 
both Z and n. The accuracy attainable remains to be determined. 

In order to define the limits within which measurement of thickness and 
refractive index can be made, observations on both continuous and skele- 
tonized cadmium arachidate films have been extended to the point at which 
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the optical distance traversed in the film by light reflected from the i 
coated chromium base is about one wavelength. This completes in roug 
outline a picture of the changes which can be brought about by the presence 
of non-absorbing films on chromium, except for the effects of birefringence, 
which have been neglected (cf. Rothen’). The cyclic nature of the 
changes can be indicated by plotting tan y exp: (iA), as in FIGURE 9. 
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Ficure 9. Effect of film refractive index on optical behavior of thick films on chromium. The data of 
FIGURE 8 are extended to include one complete optical cycle. This is represented here by plotting 
tan y-exp(i A), the origin for clean chromium being at the point O. 


It is clear from FIGURE 9 that, for each refractive index value, there will 
be a definite range of film thickness appropriate for measurement of L 
and ». Each cycle consists characteristically of a portion common to all 
films. whatever their refractive index, and an outer arc which shifts with 
changes in film refractive index. In order to be useful from the present 
point of view, experimental points must lie in the outer zone. Delimita- 
tion of the zones with a view to the setting up of experimental conditions 
suitable for the optical measurement of particular types of film is much to 
be desired but scarcely attainable at present with the calibration methods 
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used in this investigation. Diagrams of tan y exp (iA) should be available 
for several values of m, and they should be traversed by lines of equal film 
thickness to facilitate interpolation. Calculation of theoretical curves from 
the generalized Fresnel equation (Winterbottom!) has been shown by 
Mr. R. Hartman to give results which are more closely in accord with our 
experimental data for thin films than the approximate algebraic formula- 
tions considered on page 1065; the remainder of the cycle is also correctly 
reproduced. It is anticipated, therefore, that theoretical calculation will 
provide an adequate basis for the construction of calibration curves cover- 
ing a wide range of film characteristics. 

Characterization of Double or Laminated Films. The foregoing consid- 
erations indicate that the problem outlined on page 1064 can be solved 
within certain clearly defined limits by polarization optical methods. A 
single experimental point, suitably placed, is open to an interpretation in 
terms of film thickness and refractive index which is free of ambiguity, 
provided the film is known to be optically homogeneous. This is not always 
the case. Particular interest attaches to the deposition of one film upon 
another, the optical properties of which may differ considerably from those 
of the first. The ambiguity which arises when the second film is charac- 
terized merely in terms of an “optical,” or apparent, thickness has been 
considered previously! in connection with the interference method. The 
ambiguity exists equally in studies made hitherto by the polarimetric 
method." ” It is clear qualitatively from an inspection of FIGURE 9 that 
in the “outer” zone a sudden change in the refractive index of the film 
during its gradual increase of thickness will bring about a deviation from 
the loop characteristic of the material constituting the initial layers and 
that, if a sufficient amount of the second substance is added, the new 
curve will gradually approach the closed curve characteristic of films con- 
sisting only of this second substance. VaSitek!® has reduced the mathe- 
matics of the laminated film problem to that of a single film. A study by 
this method is being carried out, together with measurements on heter- 
ogeneous films. 


Summary 


The effects of thin films of proteins and other substances upon the optical 
properties of polished metals afford a valuable means of studying the 
molecular characteristics and interactions of such substances. As a pre- 
liminary to biophysical studies of this kind, we have investigated the pos- 
sibilities of the optical method by determining the behavior of multilayers 
of aliphatic acids deposited upon metal surfaces. These films provide an 
opportunity for testing the validity of the various approximate theoretical 
formulations which are available. 

Glass slides coated with chromium are examined in a polarizing assem-' 
bly suitable for the measurement of elliptically polarized light. The de- 
sired number of paired aliphatic acid monolayers is then transferred to the 
slide and the polarimetric measurement repeated. 

The relationship between the parameters of the elliptically polarized 


1080 Annals New York Academy of Sciences 


reflected light and film thickness has been determined experimentally for 
barium stearate and cadmium arachidate multilayers (refractive index 1.5) 
and for skeletonized cadmium arachidate (refractive index 1.2-1.3). The 
following conclusions have been drawn: oe 

A. The Drude approximate equations represent adequately the limiting 
slopes of the curves showing the changes of the phase parameter A and the 
amplitude parameter tan y with increasing film thickness. Significant 
deviations from the Drude curves become apparent in the case of tan wy for 
films 50 AU thick and when the thickness exceeds 100 AU in the case of A. 

B. The equations given by F. A. Lucy predict correctly that these 
deviations are in quantitative agreement with our data for A. In the case 
of tan y, they represent an over-correction. 

C. The Lucy equations imply a certain linear relationship (see equa- 
tion 6). The experimental data show a similar linearity but the constants 
of the empirical equation are not those predicted by the theory. 

D. Although the measurement of A and tan y enables film thickness and 
refractive index to be determined in principle, the experimental data on 
films less than 600 AU thick show these quantities to be so insensitive to * 
changes in refractive index that reliable values are not to be expected with 
present methods. 

E. When film thickness is comparable to wavelength, A and tan y 
undergo cyclic changes. In a certain portion of each cycle, the values 
become very sensitive both to thickness and to refractive index. It is 
possible to define the circumstances under which this fact can be utilized 
in the optical study of films by extending the experimental calibration 
curves with the aid of theoretical curves derived from the generalized Fres- 
nel equations. 

Acknowledgment is due to Miss Helen Woods for preparing the chro- 
mium-plated slides used in this investigation, to Mr. Richard Hartman 
for taking part in the work in its final stages, and to Mr. John T. Quynn 
for design of parts of the polarizing assembly. 
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GONIOMETRIC ANALYSIS OF CRYSTALS 


By J. D. H. Donnay 
The Johns Hopkins University, Baltimore, Maryland 


The methods of goniometric analysis can cope with problems of both 
determination and identification. The latter consists in checking a hy- 
pothesis as to the nature of a crystal (question: “Is X really A?”), the former, 
in finding what it is without any preconceived idea as to its composition, 
(question: “What is X?”). In either case, the first step is to obtain deter- 
minative data in the form of measured angles. 

Several kinds of angles can be measured. The most common one has 
long been the interfacial angle, usually defined as the supplement of the 
true dihedral angle between two faces. It is equal to the angle between the 
two face-normals. In present day practice, the coordinate angles of a face 
are widely used. They are defined as follows: Imagine the crystal placed 
inside a sphere of arbitrary radius, with a vertical set of parallel edges or, 
using geographical language, in the north-south direction. Then, from the 
center of the sphere, drop a perpendicular onto each face and produce it 
until it meets the surface of the sphere. The sphere is thus covered with 
an array of face-poles, the points where the face-normals puncture the 
sphere. The geographical coordinates (longitude y and colatitude or polar 
distance p) of a face-pole are the coordinate angles of the corresponding 
face. 

In addition to faces, edges also yield angles. Since a zone is defined as 
the assemblage of faces parallel to a given edge, the (plane) angle between 
two edges is an interzonal angle. The inter-edge angle between two edges 
lying in a face is usually measured in such a way that the face is comprised 
between the sides of the angle. The interzonal angle is then defined as the 
supplement of the inter-edge angle. 

The type of angles used depends on the kind of crystals to be measured 
and the equipment at hand. 


Obtaining the Data 


For large crystals (say, one cm. across or more) with rough or non- 
reflecting faces, interfacial or interzonal angles can be measured by contact 
goniometry with an accuracy that rarely attains the half-degree. A contact 
goniometer is nothing more than a 180° protractor with a moving arm 
pivoting about the center. The edge of the protractor and the edge of the 


arm are applied to the angle to be measured. The method is crude and is | 


resorted to only when other methods fail. 

If the crystal faces reflect light properly, a crystal, large or small, can be 
measured on an optical goniometer. Angles are measured with an accuracy 
of a few minutes, or even a half-minute in exceptional cases. The one- 
circle goniometer yields interfacial angles, requiring a different setting of 
the crystal for each zone to be measured. ‘The two-circle goniometer gives 
the coordinate angles of each face, one setting of the crystal usually being 
sufficient. This advantage is such that the two-circle instrument has 
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generally supplanted the one-circle model. The techniques of crystal goni- 
ometry, measurements and calculations, are readily available from a number 
of sources.!-6 

Until recently, the chief obstacle to the use of goniometric methods of 
crystal analysis had been that goniometers were not commonly available, 
particularly in this country. This difficulty no longer exists, inasmuch as 
several firms, since 1946, have built and marketed new two-circle goniometers 
(L. C. Eichner, Bloomfield, N. J.; Laboratory Associates, Belmont, Mass.’; 
Unicam, Cambridge, England). An interesting improvement of the instru- 
ment has been proposed by Terpstra.’ He places the optical axes of the 
collimator and telescope in a vertical plane (perpendicular to the horizontal 
circle), although these two optical axes were horizontal in previous designs. 
Since the arm that carries the vertical circle swings around the (vertical) 
axis of the horizontal circle, the Terpstra set-up allows a considerably 
greater range to the swinging motion. 

For small or very small crystals, and for fragments of crystals, with few 
faces or no face at all, X-ray goniometry, by means of the Weissenberg 
instrument, for instance, is indicated. This method yields absolute length 
measurements in addition to angular data. The equipment has become 

_ popular, owing to the introduction of the equi-inclination model, designed 
by M. J. Buerger® and manufactured by several firms (Otto von der Heyde, 
Cambridge, Mass.; Charles Supper, Newton, Mass.; Unicam, Cambridge, 
England). The crystal carrier of the two-circle optical goniometer should 
be interchangeable with that of the Weissenberg X-ray goniometer, and 
this specification is met by the manufacturers mentioned here. 

Finally, very small crystals (say, smaller than 0.1 mm.) can be investi- 
gated by the methods of microscope goniometry,!® in which angles are 
measured by means of the rotating stage of the microscope. Angles be- 
tween edges (interzonal angles) are measured, either in true value, if the 
edges lie parallel to the stage of the microscope, or as projected angles, if 
at least one of the two edges including the angle is oblique to the stage. 
Projected, as well as true, interzonal angles can be used to solve the crystal- 
lographic problem. No special equipment is needed, beyond a microscope 
with a graduated rotating stage. The accuracy of angular measurements 
is low, one-third or one-half of a degree. Nevertheless, microscope goniom- 
etry is a useful technique that should not be overlooked when a non-des- 
tructive method of analysis is a prerequisite. 


Using the Data for Determination 


It is obvious, of course, that goniometric methods are not designed for 
the quantitative analysis of mixtures. They can identify only pure chem- 
ical compounds occurring as crystals. 

Determinative tables in which unknown crystals can be looked up must 
be used in conjunction with goniometric analysis. The chief problem in 
the compilation of such tables lies in the choice of the determinative criteria. 
Either measured values or values calculated from measurements can be 


used. 
The Boldyrev method" aims at avoiding all calculations. The crystal 
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angles are measured and looked up as such in an index. No special orien- 
tation needs to be figured out. As Haan” puts it, the Boldyrev determina- 
tive method is crystallographic only in its aim; in essence, It 1s “purely 
administrative.’ The Boldyrev index had not been completed at the time 
of Boldyrev’s death (1946). Only part of it, the section dealing with the 
tetragonal crystal system, has been published.” ie 

Boldyrev’s efforts were an attempt at simplifying the determinative 
method of his master, Fedorov.4 The celebrated Fedorov index, Das 
Krystallreich, was completed in 1920. The determinative procedure rests 
on a proper choice of crystal setting. In Fedorov’s mind, the “‘correct.« 
setting should have structural significance. Essentially, the measured 
angles (the interfacial angles, for instance) are used to determine the shape 
of the unit cell. In the most general case, that of a triclinic crystal, five 
independent angles are necessary and sufficient to calculate the axial ele- 
ments, namely, the angles between the cell edges (interaxial angles a, 8, y) 
and the ratios of the lengths of the cell edges to one another (axial ratios 
a:b and c:b). The choice of the axes of coordinates, parallel to edges or 
possible edges of the crystal, and the choice of the parametral face (or unit 
face) intersecting all three axes, however, must be governed by a set of rules 
in order that the unit cell be uniquely determined. It was Fedorov’s aim 
to determine a unit cell that would be similar in shape to the structural 
unit cell. His rules, based on the classical law of Bravais, were inadequate 
to solve his structural problem. In so far as uniquely defining a unit cell, 
however, they were satisfactory, though complicated.® 

Another, more successful, attempt at simplifying the procedure of crystal 
determination from angular measurements is that of Barker.1® In the Barker 
method, axial elements need not be calculated; angular values are used as 
determinative criteria. A unique setting of the crystal, however, even 
though devoid of structural significance, is a prerequisite. It is obtained by 
applying a number of conventional rules designed for the purpose. Angles 
between faces with simplest Miller indices (i.e., 0’s and 1’s) are measured 
if possible, or, if not, are calculated from other measured angles. The 
determinative angle in the orthorhombic system, for instance, is the angle 
a(100):m(110), which is equal to arctan (a/b). It is convenient to be able 
to look up the angle, rather than its tangent, in the determinative table. 
The details of the method and references to the literature have been given 
elsewhere”; they can, therefore, be omitted at this point. 

After Barker’s death (1930), the compilation of his.index was continued 
at Oxford University under the editorship of R. C. Spiller and Mary W. 
Porter, with the collaboration of a number of crystallographers. The first 
volume of the index covers all crystal systems, other than monoclinic or 
triclinic, and it should appear in the very near future. It will contain 
approximately 3,300 substances. An introduction (by M. H. Hey) will 
give a full description of the Barker method and the use of the tables. Work 
on the second volume of the Barker Index is being actively carried on. 

Whereas the classical law of Bravais and the Fedorov principles that 
stemmed from it were not powerful enough to determine the structural unit 
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cell, present-day methods of X-ray diffraction can solve this problem in a 
straightforward way, and the method of morphological analysis, 1% 20 
based on the generalized law of Bravais,™” although not foolproof, leads 
to the structural unit cell (and space group) in the vast majority of cases. 
Moreover, the conventional rules governing the orientation of the unit cell 
have been revised,” ** so that it is now possible to use the structural unit 
cell for determinative purposes. 

With the help of several co-workers, I have for some years been compiling 
goniometric data in the form of axial ratios and interaxial angles. The 
determinative property is the axial ratio a/b, in the conventional setting. 
In the index now in preparation, other axial elements are also listed, together 
with absolute lengths of the cell edges (whenever they are known from 
X-ray results), space group, and specific gravity. When complete, such 
an index will make it possible to determine a crystal from either relative or 
absolute unit-cell dimensions. 

No matter which index is used, its réle in the determinative procedure is 
simply to restrict the number of possible hypotheses. In most cases, the 
determination of a crystal must be checked by means of an additional test, 
chemical, optical, etc. A combination of several methods is often fruitful. 
Optical goniometry combined with an X-ray powder pattern, for example, 

will yield the absolute dimensions of the unit cell. The latter, combined 
with a specific gravity determination, gives the weight of the cell, from 
which the molecular weight or a simple multiple of it can be calculated. 


Conclusion 


As compared with other methods of analysis, crystal goniometry occupies 
a place of its own. It cannot be expected to be universally applicable, but 
it will succeed in special cases in which other methods would fail. An 
example of such a case is that of a pseudomorphous crystal, resulting from 
the alteration or replacement of a pre-existing crystal of a different chemical 
composition, whose external form is preserved. The destroyed crystal can 
still be identified. Crystals containing inclusions (as mechanical admix- 
tures) can be determined by goniometry, since these inclusions do not affect 
the crystal form. Other advantages of goniometric analysis are fairly obvi- 
ous: the method is non-destructive; only a very small amount of material 
is required (a crystal, a fraction of one mm. in its longest dimension will 
suffice); the determinative data can be obtained in a number of different 
ways, so that various kinds of equipment can be used; the method is usually - 
fairly rapid (from about fifteen minutes to a few hours); and its specificity 
is very satisfactory, except in cases of solid solution (mixed crystals of 
similar compositions have similar angles) and homeomorphism (crystals of 
different compositions but related crystal structures may show some simi- 
larities in their angles). 

Goniometric analysis can compete with the X-ray. powder method of 
determination in dealing with crystals of organic compounds; which are apt 
to have large unit cells and low symmetry, and therefore give powder pat- 
terns characterized by a great many closely-spaced diffraction lines. 
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One limitation of goniometric analysis is its inability to cope with crystals 
of the isometric (cubic) crystal system, in which the angle between two 
faces with given Miller symbols is the same, regardless of chemical com- 
position. Even in such cases, however, goniometry is useful, as it enables 
one to identify the crystal forms (and possibly the space group), which in 
themselves are of determinative value. 
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QUANTITATIVE ANALYSIS BY MEANS OF CATALYTIC 
HYDROGENATION REACTIONS* 


By E. C. Dunlop 


Experimental Station, E. I. du Pont de Nemours and Company, Inc., 
Wilmington, Delaware 


The analytical chemist of today makes use of many different properties 
of a material in order to characterize it as completely as possible. For 
example, some chemicals, such as unsaturated organic compounds, will 
react with gaseous hydrogen, under specified conditions, with a resulting 
loss in gas volume which can be measured. Such measurement permits 
quantitative analysis by means of catalytic hydrogenation reactions. Just 
as is the case with many of the techniques used by the analyst today, cata- 
lytic hydrogenation data require skillful interpretation and are sometimes 
useful only when interpreted along with other data obtained by separate 
techniques. 

It has been known for a long time that unsaturated organic molecules 
could be hydrogenated in the presence of a catalyst and under certain 
specified conditions. This property, however, has been used only recently 
as an analytical tool. A number of papers have been published describing 
the apparatus and techniques used in carrying out the reaction. To men- 

tion some of these, articles by Prater and Haagen-Smit,! Johns and Seiferle,? 
Shriner,® Hallet,* Noller and Barusch,® and Joshel,® have appeared. 


Apparatus 


We have studied these various types of apparatus and have constructed 
similar but modified equipment. FicurRE 1 shows a diagram of our present 
equipment. This equipment is more closely related to that described by 
Lloyd M. Joshel® than to any of the others. 

Interesting points about the equipment are listed below. The primary 
modifications which we have made are described in sections three and five. 

1. The stopcocks are all oblique bore and Burrell’s #2 heavy stopcock 
grease is used. 

2. The reaction flask has a male grind and the buret system the female 
grind, thus helping to avoid contamination of the reaction flask with stop- 
cock grease. 

3. The reaction flask is immersed in a constant temperature bath and 
the burets are jacketed in order to circulate constant temperature liquids. 

4. The connecting tubing is capillary tubing with dimensions of 8 mm. 
O. D., 3 mm. I. D. 

5. The reaction flask has a number of sealed side arms for holding the 
sample or catalyst so that it can be added to the solvent as needed. The 
sample or catalyst is placed in a small glass boat which has a piece of iron 
sealed in one end. The boat is pushed into the flask when desired by means 
of a solenoid. Additional portions of the sample or catalyst can be added 
by having several boats in the different side arms. 


*Contribution No. 265 from the Chemical Department, Experimental Station, E. I. du Pont de Ne- 
mours and Company, Inc., Wilmington, Del. 
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Procedure 


We have standardized our procedure and technique, as have others, with 
the use of maleic acid. A pure sample of maleic acid having a melting 
point of 130.5°C. was used. This unsaturated dibasic carboxylic acid reacts 
with hydrogen to give succinic acid, one mole of Hy per mole of maleic acid 
being absorbed. In confirmation of previous results reported in the litera- 
ture, it was found that a sample size of approximately 100 milligrams and 
ten per cent that amount of catalyst gave the best operating quantities 
for our apparatus. 


‘CATALYTIC HYDROGENATION APPARATUS 


I PUYUTOVUONIVOVETOVETO ON ON ON OV OO EY ITED | ~) 


fa PURTRNOROTOTEVETOT TOT NTV TED 


4 -~~-REACTION FLASK 0 ~-- MANOMETER & SCALE G& OSMERGSs INLET 


B----- CONSTANT TEMP. E & E- CONSTANT TEMP. hats 
LIQUID BATH LIQUID INLET 
& C- BURETTES F &F- CONSTANT TEMP. S---- MAGNETIC 
LIQUID EXIT STIRRER 
FIcureE 1. 


The procedure adopted for carrying out analyses may be outlined as 
follows: 

The apparatus should be clean and dry, and sufficient mercury should 
be in the burets and leveling bulbs. 

The circulating pump should be turned on to circulate the constant 
temperature liquid through the buret jackets. 

The reaction flask should be clean and dry and a small amount of grease 
applied to the male grind. Care should be taken to keep the grease from 
the top of the grind. Ten milligrams of catalyst and 15 milliliters of solvent 
are placed in the bottom of the flask and a glass enclosed magnet is also 
inserted. Approximately 100 milligrams of the sample is accurately weighed 
into the sample boat and placed in one side arm of the flask. The flask is 


then attached to the apparatus and the grind well seated and attached 
with spring clamps. 
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The constant temperature bath is then raised to surround the reaction 
flask and the magnetic stirrer is turned on. The system is then alternately 
evacuated and flushed with hydrogen four or five times and finally left 
filled with hydrogen. The leveling bulbs are adjusted so that the burets 
read 0.0 milliliter at atmospheric pressure. This is accomplished by venting 
or adding hydrogen. After taking the initial reading, the boat containing 
the sample is forced into the reaction flask by the solenoid and a positive 
pressure of 15-20 centimeters of mercury is applied by raising the leveling 
bulbs. The barometric pressure is recorded at the beginning and at the 
end of each determination. Hydrogenation is continued until no further 
decrease in pressure is noted on the manometer. This usually requires 
about 1 to 15 hours. When hydrogenation is apparently complete, the 
leveling bulbs are adjusted to give atmospheric pressure in the apparatus 
and the buret reading is taken. The leveling bulbs are again raised to 
apply pressure of about 15-20 centimeters of mercury on the system for 
about 10 minutes. The system is then brought to atmospheric pressure and 
the burets read a second time. If this reading checks the previous one, the 
reaction is complete and the volume of hydrogen used is converted to stand- 
ard conditions and the unsaturation value of the sample is calculated. 

We have adopted the practice of reporting data in terms of grams of 
hydrogen per gram of sample. The data are used in several ways and this 
appears to be the most useful single form for all concerned. We have used 
the following catalysts: PtO2, Pd-on-charcoal, nickel catalyst W-6 (Adkins- 
Billica, 1948), and Baker Colloid 46 (Palladium). We have used the follow- 
ing solvents: methanol, ethanol (95 per cent), ethanol (abs.), water, glacial 
acetic acid, ethyl acetate, and dioxane. 


Results 


We will now take up the results which we have obtained. 

First, with maleic acid, the following table shows the results obtained 
with various sizes of samples, different catalysts and different solvents. 

TABLE 1 indicates that 95 per cent ethanol, acetic acid, and dioxane 
solvents give the best results. Also, it may be seen from the table that 
either platinum oxide or palladium-on-charcoal catalyst can be used satis- 
factorily. 

This has not been an exhaustive study, and it is not intended by these 
data to prove that the other solvents such as ethyl acetate and water are of 
no value. The purpose of this experiment has been to find the most satis- 
factory experimental conditions under which to carry out our analyses. 

We have used this apparatus to determine the amount of unsaturation 
in high vapor pressure samples such as diisobutylene. In this experiment 
we used ice water in the circulation system, and the results, which check 
well with the theoretical value, are shown in TABLE 2. 

Some work has been done with unsaturated ketones, and it is possible to 
determine the carbon-carbon unsaturation using platinum oxide or palladium 
on-charcoal catalyst in 95 per cent ethanol without reducing the carbonyl 
group. The carbonyl group can be reduced with nickel catalyst, prepared 
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TABLE 1 
HypROGENATION OF Mateic Acip 
Theory—Grams H2/g. Sample = 0.017 


H2/g. 
Sample Wt.-mg. Catalyst-mg. Solvent wee cet 
100.8 10 PtOz Ethanol (95%) 0.017 
96.6 10 PtOs Acetic acid 0.017 
84.7 10 PtO2 Dioxane 0.017 
105.4 10 PtOz Methanol 0.016 
90.8 10 PtOz Ethanol (abs.) 0.018 
87.6 10 PtOz Ethyl acetate 0.015 
85.6 10 PtO, Water 0.013 
88.4 10 Pd on Charcoal Acetic acid 0.017 
85.0 10 Pd on Charcoal Dioxane 0.017 
85.3 10 Pd on Charcoal Ethanol (95%) 0.017 
10 Baker Colloid 46 Ethanol (95%) 0.016 
TABLE 2 
HYDROGENATION OF DIISOBUTYLENE 
Theory—Grams H»/g. Sample = 0.0179 
Sample wt.-mg. PtO2 Catalyst wt.-mg. Grams H2/g. sample 
102.0 17 0.0170 
86.0 10.0 0.0178 
113.6 9.1 0.0168 
TABLE 3 
HYDROGENATION OF UNSATURATED KETONES 
Grams H2/g. sample 
Sample Catalyst Solvent ui : 
Found Theory 
Methyl vinyl ke- Pd-on-charcoal | 95% Ethanol | 0.0320 | 0.0288 (1C=C) 
tone (A) A = 0.0310 
Acetone Ni (W-6) 95% Ethanol | 0.0338 | 0.0348 (1C=0) 
; 0.0340 
cone nel ke- PtO, 95% Ethanol | 0.0240 | 0.0288 (1C—=C) 
one (B) r 0.0232 
Methy vinlyl ke- Ni (W-6) 95% Ethanol | 0.0455 | 0.0576 (1C=C) 
tone (B) 0.0454 (1C=0) 
TABLE 4 


HYDROGENATION OF STYRENE 


Grams H2/g. sample 


Sample Catalyst Solvent 
Found Theory 
aoe a ee 61°C at | PtO. | 95% Ethanol 0.0186 | 0.0194 1(C=C) 
é 0.0176 
N%—1.5439 PtO: | 95% Ethanol + 4] 0.0751 | 0.0775 4(C=C) 
ml. HCl 0.0694 


95% Ethanol + 1} 0.0752 
ml. HCl 


6eR+F7+ ooo 
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by the procedure described by H. Adkins and H. R. Billica,? in 95 per cent 
ethanol, which is made slightly alkaline. 

TABLE 3 shows data for the reduction of the C=C bond in methyl vinyl 
ketone, the reduction of the carbonyl group in acetone, and the reduction 
of both the C=C bond and C=O bond in methy] vinyl ketone. 

We obtained results that are close to the theoretical values for the first 
two samples above. The results for the methyl vinyl ketone sample (B) 
using PtO2 catalyst and calculating on the basis of reduction of'one C=C 
bond indicate 81 per cent MVK and using Ni catalyst (W-6) and calculating 
on the basis of reduction of both the C=C and C=O bonds indicate 78 
per cent MVK. 

The next experiment under discussion will be the analysis of a sample 
of cyclooctatetraene. This is an example of how more than one technique 
is involved in chemical analysis and also shows how intelligent interpreta- 
tion is required from the analyst to both obtain and report the data. Ultra- 
violet absorption data were used in conjunction with the hydrogenation 
results. 

The sample reacted with hydrogen in the presence of PtO, catalyst in 
95 per cent ethanol to absorb 0.075 and 0.076 gram H»/gram of sample. 
This corresponds to 97.5 per cent cyclooctatetraene calculated on the basis 
of four C=C bonds. Ultraviolet absorption data on the original sample 
indicated the presence of styrene. The absorption spectrum of the reduced 
sample was similar to that of ethylbenzene and, calculating the concentra- 
tion of ethylbenzene, based on the absorption coefficient as given by the 
API (American Petroleum Institute) data for ethylbenzene, we got a value 
of approximately 5 per cent for the styrene in the original cyclooctatetraene 
sample. 

A known sample of styrene was hydrogenated and the results in TABLE 4 
show that the above experiment was correctly interpreted in that only the 
ethylenic bond was hydrogenated under the conditions used. Also included 
in TABLE 4 are data for the hydrogenation of the aromatic portion of the 
styrene. By the addition of acid, we were able to completely hydrogenate 
the sample of styrene. 


Summary 


A reaction flask has been described which has sealed side arms used for 
holding the sample until it is needed. This enables the operation to be 
carried out in a closed system, eliminating one leak hazard and simplifying 
the operation. Data have been shown for maleic acid in various solvent 
and catalyst combinations, and it is indicated that 95 per cent ethanol, 
acetic acid, and dioxane are excellent solvents with either PtO, or Pd-on- 
charcoal catalyst. 

Unsaturated ketones (methyl vinyl ketone) have been analyzed by selec- 
tively hydrogenating the C=C and the C=O, using PtO: or Pd-on-charcoal 
catalyst for the C=C and Ni W-6 catalyst for the C—O. 

The analysis of a sample of cyclooctatetraene containing some styrene 
has been made in conjunction with ultraviolet absorption spectrum analysis. 
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Cyclooctatetraene and the ethylenic group of styrene are hydrogenated in 
neutral solution with PtO2 catalyst. The aromatic portion of the sample 
is hydrogenated along with the ethylenic bond when carried out in acidic 
alcoholic solution using PtO, catalyst. 
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ANALYSIS OF MIXTURES BASED ON RATES OF REACTION 


By T. S. Lee* and I. M. Kolthoff 


School of Chemisiry, University of Minnesota, Minneapolis, Minnesota 


The analysis, by conventional means, of mixtures of two constituents 
that have similar chemical properties is often difficult or impossible. Mix- 
tures of two organic compounds that contain the same functional group 
(e.g., two esters or two aldehydes) are frequently encountered and present 
special difficulties to the analyst. For example, mixtures of homologous 
organic compounds or of certain ortho- and para-isomers usually cannot 
be analyzed by classical chemical methods. Methods that depend on 
physical properties of the compounds are usually not suitable either. In 
illustration of this point, ultraviolet spectrophotometry cannot be used to 
analyze a mixture of homologous ketones because the spectrum of each 
ketone is determined chiefly by the carbonyl group and is affected only 
slightly by the nature of the alkyl group substituents. Infrared spectro- 
photometry, which is capable of distinguishing between homologous or- 
ganic compounds, requires expensive equipment and is frequently not 
applicable if the sample is available only as a dilute solution in water or 
in certain other common solvents. (On the other hand, infrared analysis 
is ideal if the apparatus is available and if the sample is available in a 
suitable form.) 

It frequently happens that such mixtures of organic compounds can be 
analyzed by a method based on differences in rates of reactions of the com- 
ponents of the mixture with a given reagent. The effect of a substituent 
on the rate of reaction of a reactive group with a particular reagent is 
usually pronounced. Propylene, for example, reacts with perbenzoic acid 


‘at only one-thirtieth the rate at which 2-butylene reacts (see below). 


Methyl alcohol reacts with hydrobromic acid nearly six times as rapidly 
as ethyl alcohol does (solvent phenol, 80°). It should be pointed out, 
however, that, if the substituents of the functional group are quite similar, 
as for example, -butyl and m-amyl, the method of analysis based on reac- 
tion rate cannot be applied. Furthermore, the method is generally limited 
to mixtures of two components. 

Highly specific methods for the analysis of mixtures based on rates of 
reaction have been described previously in the literature.~’ Most of these 
methods were developed for the investigation of specific physical chemical 
problems and are not convenient as general analytical procedures. They 
differ greatly from the proposed general method. In order to understand 
the principles underlying the proposed method of analysis, consider two 
compounds, A and B, that react with a reagent R at different rates. Sup- 
pose that the reaction of A with R is faster than the reaction of B with R, 
at a given constant temperature. Mixtures of A and B can be analyzed 
by the following procedure. A reaction solution is prepared in which the 
the total molar concentration of A + B is specified and in which the con- 


* Present address: Department of Chemistry, University of Chicago, Chicago, Illinois, 
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centration of R is also specified. The reaction is allowed to proceed for 
a specified time at a given constant temperature and then the amount of 
A + B reacted is determined. It is apparent that, if the mixture is com- 
posed chiefly of A, the more reactive compound, a relatively large fraction 
of A + B will have reacted after the specified reaction period. On the 
other hand, if the mixture is composed chiefly of B, the less reactive com- 
pound, a relatively small fraction of the mixture will have reacted. It 
should be noted that the method is an indirect one and possesses the dis- 
advantages common to other indirect methods of analysis. The deter- 
mination, therefore, of a small amount of one component in the presence 
of a large amount of the other is generally difficult. 

A concrete illustration of the method is the procedure for the analysis of 
mixtures of ethyl and isopropyl acetates. (This procedure is based on 
the rates of reaction of the two esters with hydroxyl ion.) First, the total 
ester content of the sample is determined. Second, a reaction solution is 
prepared at 0°C in which the total concentration of ester is 0.0500 M and 
the initial concentration of alkali is also 0.0500 M. The mixture is allowed 
to stand for 33 minutes at 0°C and the concentration of hydroxyl ion re- 
maining is then determined by acidimetric titration. Third, the concen- 
tration of hydroxyl ion remaining after the reaction period of 33 minutes 
is compared with a calibration curve that gives the initial composition of 
the mixture directly. Thus, two titrations (the determination of total 
ester content of the sample and the determination of the hydroxyl ion 
remaining at the end of the reaction period) are sufficient to establish the 
original composition of the mixture. 

In general, the reactions of the two components of the mixture, A and B, 
with the reagent, R, may be called competitive reactions and be repre- 
sented as follows: 


A + R-— products (1) 
B+ R=— products (2) 


So far as the analysis of mixtures is concerned, the most interesting and 
useful type of competitive reactions is that in which A and B are organic com- 
pounds that contain the same functional group. In such a case, the two 
reactions, (1) and (2), generally proceed by identical mechanisms and differ 
only in specific rates of reaction. Such a pair of reactions we propose to 
call homo-competitive reactions, the prefix homo- referring to the fact that 
both compound A and compound B contain the same functional group. 


The term homo-competitive reactions is extended to include first order 
reactions of the type: 


A — products (3) 
B — products (4) 


in which A and B are organic compounds containing the same functional 
group. | An example of first order homo-competitive reactions is the de- 
composition of formaldehyde bisulfite and acetaldehyde bisulfite. 


eat 
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In addition to the homo-competitive reactions mentioned above, there 
are many other types of competitive reactions. For example, the simul- 
taneous reactions of bromine with an ethylenic bond (addition) and with 
a methylene group (substitution) might be called hetero-competitive re- 
actions, since the reagent (bromine) is reacting with two different functional 
groups. Hetero-competitive reactions are not considered in this paper. 
They belong to the category in which one or more “‘side-reactions” occur, 
in addition to the main reaction with the reagent. 


Mathematical Treatment 


Although mathematical investigations of certain aspects of competitive 
reactions have been reported in the literature,®: *: !° these investigations are 
of no use for the present purposes and will not be considered further. 

It is possible to treat the kinetics of two homo-competitive reactions in 
an exact theoretical way. On the basis of this treatment, it is possible 
to calculate calibration curves for the analysis of mixtures from the values 
of the specific rate constants of the two homo-competitive reactions. (Al- 
ternatively, the calibration curve may be established empirically by analy- 
zing mixtures of known composition.) 

In order to calculate a calibration curve for the analysis of mixtures of 
A and B, it is necessary to use an equation relating the concentrations of 
A and B remaining to the reaction period. Equations will be derived for 
first order homo-competitive reactions and also for second order homo- 
competitive reactions. 

First Order Reactions. ‘The concentration of A remaining at any time, #, 
is given by the expression: 


a=ae'? (5) 


In this equation, a@ is the concentration of A at any time, ¢; dp is the initial 
concentration of A; and k’ is the specific rate constant for the reaction of A. 
A similar expression gives the concentration of B at any time, ¢. The sum 
of concentrations of A and B is, therefore, given by the expression: 


Gb = He + he * (6) 


It is apparent from equation (6) that if the values of a + b, a + bo, #, F’, 
and R” are known, the individual values of a and 6b can be found (see below). 

Second Order Reactions. If A and B react with a reagent by second order 
reactions (see reactions 1 and 2), the rates of disappearance of A and B 
will be given by the following expressions: 


2 ay aap ae (7) 
Breet kh’ br (8) 
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In these equations, k’ and k’’ are the specific rate constants for the reactions 
of R with A and with B, respectively, and r is the concentration of the reagent | 
R at any time, ¢. Equations (7) and (8) can be combined to give: 


db kb 
ag eS. 9 
da ae (9) 
Integration yields: 
b Rk” a 
So peste = 10 
log 5 = G log = (10) 


or 


b - ketlke (: y 
ia |e = {— 11 
bo (<) ao ( ) 


in which a represents the ratio of rate constants, k’’:k’. (In the following 
derivations, it will be convenient to assume that ’ is greater than k’’; that 
is, the reaction of R with A is faster than the reaction of R with B. The 
value of a will therefore always be less than unity.) It is seen from equation 
(11) that the fraction of B remaining in the reaction mixture at any time 
is equal to the fraction of A remaining, raised to the power a. Calibration 
curves for the analysis of mixtures cannot be calculated directly from equa- 
tion (11) because time, ¢, does not appear in it. 

The derivation of an equation in which the concentration of A or B or R 
is expressed as a function of time is not simple. Equation (7) can be re- 
written 


—— = k' al(m — am — hh) tat J] (12) 


bs da 
a[(ro AY ae bo) +a-t+ b] 
The expression for 5 given in equation (11) is now substituted in (13). 


da 
bo 
a (r — a — bh) tat ao" | 
ao 


= k' dt (13) 


See (14) 


It is necessary now to integrate equation (14). Because of the complexity 
of the integration, only the special case, ro = ao + bo, will be considered. 
Equation (14) can also be integrated for the general case, ro ~ ap + by. 
The integrated equation for the general case is so complex, however, that 
it is of little practical value. 

If 79 = a9 + bo , equation (13) becomes: 


da P 
Cr cS 
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The left side of equation (15) is a function of a alone since: 


tae oe (16) 


(See equation 11.) Integration of both sides of equation (16) yields: 
da 


p= [ “i+3) (17) 


The substitution, y = a is made. In view of equation (16): 


Therefore: 


a\ 1/(a—1) 
ao 1 (2—a) /(e—1) 
Gh, | == d 
. (#) (@ = 1)” ‘ 


A at 2/(a—1) ae 
= | — 
bo 


Substituting these values in equation (17) gives: 


Y (by \VD 1 (yal 2) 
Ki = | & ——— 18 
yo \@o a—l 1+ y y ( ) 
b 1/(a—1) 1 
We let Cw @ | ( ) (19) 
ao a—i1 
a 


It has been assumed that a is always less than one (see above). The value 
of B is also less than one if the value of a is less than}. Equation (18) 
becomes 

y yh d 

y0 i+ My Z 


Another substitution, « = 1 + y, is made: 


oe See i vi. 
x0 x 


—kt=C (21) 


-cf'- Ae, ae (22) 
Pie nee sf ’ = =" de 
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B . . 
The term (1 — *) can be expanded by the binomial expansion: 
x 
ae fe p-2 , BIB—1) ps 
—R't = ce E = 8x e- ee x dx 


- a a Ree 8 al 
-Wi=c|5-ghoe Tone 2) by 


bo es 
inwhich B= ¢-"—33 @ = #"/R and # = Lty=t+0o=14(%2)q 


and C is defined by equation (19). 

Equation (23) is of considerable practical importance in the analysis of 
mixtures. Calculations made from this equation are given below. 

It should be mentioned that, if the initial concentration of R is much 
greater than the sum of the initial concentrations of A and B (ie., if ro > 
ao + bo), the reactions of A and B will be pseudo-unimolecular. In this 
case, the concentrations of A and B remaining after a given reaction period 
are represented by the following equation: 


at P= net i 4 bye cs (24) 


The derivation of equation (24) is similar to that of equation (6). 

Optimum Reaction Period. The calibration curve for the analysis of 
mixtures is a plot of the percentage of A + B reacted after a given reaction 
period versus the initial composition of the mixture. The slope of the cali- 
bration curve, and therefore, the accuracy of the analysis, are dependent on 
the choice of the reaction period. For every pair of homo-competitive re- 
actions, there exists an “optimum reaction period.” The calibration curve 
corresponding to the optimum reaction period has a greater slope than the 
calibration curve corresponding to any other reaction period. Conse- 
quently, the maximum accuracy in the analysis of reaction mixtures is ob- 
tained by employing the optimum reaction period. For a pair of first order 
homo-competitive reactions, the optimum reaction period is calculated from 
the equation: 


In (R'/k” 
lopt. = ee, (25) 


For a pair of second order homo-competitive reactions, if 79 = do + bo: 


lopt: = ee (26) 


For a pair of second order homo-competitive reactions, if 75 > ay + bo: 


In (k'/k!’) 
roll = oe 


lopt. = 


(The aaa of equations 25-27 will be given in a subsequent publi- 
cation. 
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Numerical Calculations 


Calibration curves for the analysis of mixtures of A and B can be found 
empirically by analyzing mixtures of known composition. Alternatively, 
the calibration curve may be calculated from equations (6), (23), and (24), 
provided that the values of the specific rate constants of the two reactions 
are known. The calculations are sometimes tedious but can be greatly 
shortened if the following procedures of calculation are used. 

First Order Reactions. The calibration curve for a pair of first order homo- 


competitive reactions is found by plotting values of . - ; (calculated from 
0 0 
equation 6) versus values of . r b Inasmuch as k’, k’’, t, ro , and (ay + Bo) 
0 0- 


are constant in the analysis of mixtures, it is apparent that the calibration 
curve for first order reactions is a straight line. 


TABLE 1 


CALCULATION OF CALIBRATION CURVES FOR First ORDER 
HOoMO0-COMPETITIVE REACTIONS 


a (% (A + B) Remaining) 10% a\(% (A +B) Remaining)i0% 
t 19 57 

B 13%5 67.5 

1 8 78 

ar Sad 87.5 


If the optimum reaction period (see above) is employed, a simplified 
method of calculating the calibration curve can be used. Combination of 
equation (25) with equation (6) results in the following equations: 


= = (% of (A + B) remaining) 10% 4 = ee (28) 
0 

b . . a ‘rat | 

iad (% of (A + B) remaining) sop = pa (29) 
0 


in which the subscript 10074 refers to a “mixture” with composition 100% A, 
and a is the ratio of specific rate constants k’’/k’. Thus, if the optimum 
reaction period is employed, the calibration curve for first order reactions 
can be found by joining with a straight line the points calculated for 100% A 
(equation 28) and for 100% B (equation 29). The results of the calcula- 
tions of calibration curves for four typical pairs of first order homo-competi- 
tive reactions are given in TABLE 1 and in FIGURE 1. It should be noted 
from equations (28) and (29) that, if the optimum reaction period is em- 
ployed, the calibration curve for first order reactions is determined solely 
by the value of a, the ratio of rate constants. 

Second Order Reactions. Calibration curves for second order homo-com- 
petitive reactions are not, in general, straight lines, although usually they 
are nearly straight. The ease with which the calibration curves can be cal- 
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culated is dependent on the ratio of the two second order rate constants. 
If this ratio is greater than 25, the calibration curves can be calculated in 
a simple manner, as will be described below. If the ratio Is less than 25, 
however, the calculation is more difficult and can be performed only for two 
special cases: 79 > do + bo andro = ao + bo. | 

(a) ro > ao + bo. For this case, the reactions of A and B with R are 
pseudo-unimolecular. Consequently, equation (24) or equations (28) and 
(29) may be used in calculating the calibration curves, as described above. 
It should be noted that the calibration curves of FIGURE 1 apply in this case. 

(b) ro = ao + bo. The equation representing the rate of reaction of 


FicureE 1. Calculated calibration curves for analysis of mixtures (first order reactions). 


A and B with the reagent is equation (23) above. Equation (23) can be 
represented in the following simple form: 


—k't = C(f — fo) (30) 


In this equation, C = (2) ns = , f is the value of the expression in 
. = 

brackets in equation (23) at any time, /, and fy is the value of the same 

expression initially, at = 0. The numerical computations involved in the 

calculation of the curves from equation (30) and the values of the two rate 

constants are tedious but may be shortened if the following procedure of 

calculation is adopted: 

Step 1. The values of C (equation 30) are calculated for a series of hypo- 
thetical mixtures containing 0 per cent A, 10 per cent A, 20 per cent A, etc. 
The results of such calculations are given in TABLE 2. In this particular 
example, it is assumed that the value of a, the ratio of rate constants, is 
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0.254. (This is the ratio of rates of saponification of isopropyl and ethyl 
esters at O°C. The theoretical calibration curve calculated below is com- 


; ; TABLE 2 
CALCULATION OF VALUES OF C (a = 0.254, ro = ao + bo = 0.05 M) 


Initial composition ee TS ee aon : C 
of mixture (%A 
f (7A) ao (M) bo (M) 
0 0 0.050 0 
10 0.005 .045 —14.10 
20 010 ,040 —20.88 
30 015 .035 —28.70 
40 020 .030 —38.86 
50 025 .025 —53.62 
60 030 .020 —76.80 
70 035 .015 —119.3 
80 040, .010 —214.4 
90 .045 .005 —566.5 
100 .050 0 
TABLE 3 


CALCULATION OF VALUES OF f, EquaTIon 30 (a2 = 0.254) 


b/a ae Ist term* | 2nd term | 3rd term | 4th term | 5th term f, Sum 
of terms 
0.20 1.20 3.125 0.458 0.050 0.014 0.006 3.653 
0.25 1225 3.169 0.446 0.047 0.013 0.005 3.679 
0.30 1.30 3.211 0.434 0.044 0.012 0.004 3.706 
0.35 1235 3.253 424 O41 O11 .004 3.732 
0.40 1.40 3.294 A414 .039 .010 .003 3.759 
0.50 1.50 3.372 395 .035 .008 .003 3.812 
0.60 1.60 3.447 .379 .031 .007 neg. 3.863 
0.80 1.80 3.587 .350 .026 -005 i 3.968 
1.0 2.0 3.718 327 .021 .004 - 4.070 
iy) Bed) 3.841 307 .018 .003 4.170 
18S 2) 0) 4.012 282 .015 .002 ws 4.311 
1.75 2S 4,144 .265 013 002 4.423 
Z 3 4.270 .250 O11 001 “s 4.533 
3 & 4.708 .207 .007 001 iy 4.922 
4 5 5.084 .179 .005 neg. : 5.267 
5 6 5.404 .158 003 ss 5.566 
6 7 5.698 143 .003 re Bs 5.843 
7 8 5.962 131 .002 we t 6.095 
8 9 6.206 121 002 ve H 6.329 
10 11 6.646 .107 .001 ¢ & 6.754 
15 16 7.548 .083 neg. : # 7.631 
20 21 8.282 .069 es 8.351 
25 26 8.925 .060 “ os SS 8.985 
30 31 9.457 .054 a me Ss 9.511 


: 5 Rie ft 
* ‘4st Term” refers to the first term of the infinite series of equation (23); #.e., B af 


pared with the experimental results of the analyses of mixtures of isopropyl 
and ethyl esters in FIGURE 3.)* 


* The specific rate constant for the reaction of isopropyl acetate with hydroxyl ion is 0.303 1. mole“ min. 
at 0° and that of ethyl acetate with hydroxyl ion is 1.19 1. mole min.. Substitution of these values in 
equation (26) yields 33.3 minutes as the optimum reaction period at 0°. This reaction period is used in 
analyses of mixtures Of the two esters (see FIGURE 3). 
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Step 2. The values of f (equation 30) are calculated for several different 
ratios of b:a. The results of such calculations are given in TABLE 3. It will 
be observed that, when 0/a is less than about 0.5, five terms of the infinite 
series that appears in equation (23) must be calculated. For greater ratios 
of b/a, fewer terms are needed to yield a sufficiently accurate value of f. 


FicurE 2. Relation between f of equation 30 and b/a. (a = 0.254) 


TABLE 4 
CALcuLaTION OF b/a (a = 0.254, k’ = 1.19 min“! M~}, ro = do + bo = 0.05 M) 


Composition ap 
of mixture °) So* aot ft a 
( LA aso a 
0 a —* es — —— 
10 9.000 6.542 2.811 9.353 28.4 
20 4.000 5.267 1.898 7.165 12.15 
30 2.333 4.675 1.381 6.056 6.81 
40 1.500 4.311 1.020 555k 4.22 
50 1.000 4.070 0.739 4.809 2.665 
60 0.6667 3.898 0.516 4.414 1.730 
70 0.4296 3.774 0.332 4.106 1.073 
80 0.2500 3.679 0.185 3.864 0.600 
90 0.1111 — _ — _ 
100 — a — — — 


* Obtained from plots similar to FIGURE 2, but larger scale. 


+ The value of ¢ 1s taken to be 33.3 min., the optimum period of reaction for mixtures of ethyl and iso- 
propyl acetates. 


¢ From equation (31). 
§ From plots similar to FIcuURE 2. 


Also, the smaller the value of a, the more rapidly the series converges, and, 
consequently, the fewer the terms needed. In other words, oaly a few 
terms need be calculated if the ratio of one of the reactions is very much 
greater than the rate of the other. 


The values of f are plotted against b/a in ricURE 2. In order to obtain 
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the necessary accuracy in the steps that follow, it is necessary to plot the 


values of b/a of TABLE 3 on a very much larger scale than is shown 
in FIGURE 2, 


TABLE 5 


CALCULATION OF % or (A + B) ReMAtninG (a = 0.254, k’ = 1.19 min™! M7“, 
To = do + by = 0.05 M) 


Initial com position b % of (A + B) 
of mixture (YA) @ < remaining 
0 — = = 66.5* 
10 28.4 0.00107 0.0304 63.0 
20 12.15 .00226 .0274 59.4 
30 6.81 .00357 .0243 55.8 
40 4.22 .00500 .0211 Spe 
50 2.665 .00672 .01790 49.2 
60 1.730 .00837 .01448 45,7 
70 1.073 .01023 .01097 42.4 
80 0.600 .01240 -00744 39.7 
90 _ — = —_ 
100 — _ — lao 


* The % of (A + B) remaining in “mixtures” containing initially 0% A and 100% A was calculated simply 
by means of the ordinary second order rate law equation (33a). 


Ficure 3. Analysis of mixtures of ethyl and isopropyl acetates at 0°C. (Curve is theoretical, points are 
experimental. 


Step 3. The value of fy is found from, the plot of f vs. b/a. This is illus- 
trated in TABLE 4, column 3. The value of f is then calculated from the 
known values of fo, 2’t, and C: 

k't 
f= fo — Cc (31) 
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(Cf. equation 30). The results are given in column 5 of TABLE 4. The 
value of b/a corresponding to the value of f is then obtained from the plot 


of f vs. b/a (FIGURE 2). z. 
ca 4. The values of a and 8 are calculated from the value of b/a by 


means of the following equation: 


Paes (2) = Gre de! (32). 
bo a 


(Equation 32 is derived from equation 11). The results of such calculations 
are given in TABLE 5 and in FIGURE 3 (curve is theoretical, points are ex- 
perimental, see below). It will be noted that the calibration curve is nearly 


E) rr) C) c) 100 


° 
Anes i. Calculated calibration curve for analysis of mixtures (2 = 0.04, k’ = 1 min“! M+, ro = ag + 


a straight line. (This is also true if the difference in rate constants of the 
two reactions is large. See FIGURE 4.) 

The methods described above are useful in calculating calibration curves 
for second order homo-competitive reactions if the ratio of rate constants 
is less than 25 and if ro = ap + bo or 79 > Go + bo. If the ratio of rate con- 
stants is greater than 25, however, a more convenient method of calculation 
can be employed. In the simplified method, it is assumed that the reaction 
of A (the faster reacting component of the mixture) with the reagent R is 
represented by the ordinary equation for a second order reaction:* 


. i yas : F ; 
ei = Lp, 0 (to — +) (33) 
(70 — ao) 7% a 
& *If ro = ao, equation (33) is indeterminate. In this special case, the equation: 
aE 
k't=--—— (33a) 


, a a 
must be used instead of equation (33). : 
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Thus, the reaction of A with R is assumed to proceed as though B were ab- 
sent. Incalculating the amount cf B remaining after the specified reaction 
pericd, it is assumed that the effective initial concentration of reagent is 
ro — aoratherthan 7). That is, it is assumed that A reacts instantaneously 
with the reagent and so reduces the effective concentration of reagent from 
19 to 7) — a9. Consequently, the following equation, which is analogous to 
equation (33), can be written: 


ieee. 1 bo (ro — ao) — bo + 
tal (ro — do) — bo ig la — a) b | oP 


An illustration of the use of the simplified method of calculating calibra- 
tion curves is given in TABLE 6. It is arbitrarily assumed that k’ is 1.00 
min“ M-‘and that k”’ is 0.04 min-! M-! (therefore, ais 1/25), and that 77) = 


TABLE 6 


CALCULATION OF PER Cent oF (A + B) REMAINING BY THE APPROXIMATE METHOD 
(a = 0.04, k’ = 1 min“ M“}, ro = ao + bo = 0.05 M) 


Initial Concentrations of A and B 
' 0 Initial concentrations (M) remaining at end of reacin. % of 
om position iod (M) (A+ B) 
a ay ae wk remaining 
(% A) ao bo a b 
0 0 0.050 0 0.04167 83.3 
10 0.005 045, 0.000050 .03814 76.4 
20 .010 .040 .000145 .03448 69.3 
30 .015 .035 .000315 .03070 62.0 
40 .020 .030 .000610 .02679 | 54.8 
50 025 .025 .00107 5022735 a 47.0 
60 .030 .020 .00177 .01852 40.6 
70 .035 .015 .00278 .01415 33.9 
80 .040 .010 .00416 .00962 27.6 
90 045 .005 .00601 -00490 21.8 
100 .050 0 .00833 0 16.7 


a + b) = 0.05 M. The optimum reaction period for the analysis is com- 
puted from equation (26), and found to be 100 min. In column 4 of TABLE 
6 are given the results of the calculation of concentration of A remaining 
(equation 33) and in column 5 are given the results of calculation of the 
concentration of B remaining (equation 34). The results of calculation of 
the % of (A + B) remaining (column 6) are accurate to at least a few tenths 
of a per cent. The calibration curve is shown in FIGURE 4. 


Results of Analyses of Mixtures 


It has been found that the homo-competitive method of analysis is prac- 
tical if the rates of the two homo-competitive reactions (either first order 
reactions or second order reactions) differ by a factor of four or more. The 
following applications of the method have been made: 

Analysis of Mixtures of Esters. The reaction of esters with hydroxyl 
ion (saponification) is a second order reaction. ‘The relative rates of saponi- 
fication! 2: 18 of methyl, ethyl, isopropyl, and tertiary butyl acetates at 
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25° are in the ratio 119:72:17:1. The relative rates of saponification of 
methyl formate, acetate, monochloracetate, and dichloracetate at 25° are 
in the ratio 223:1:761:16,000. From these values, it is seen that the analy- 
sis of many pairs of esters is possible. As an example of analysis of com- 
pounds that do not react with the reagent at greatly different rates, mixtures 
of ethyl acetate and isopropyl acetate were analyzed and the results were 
compared with the calibration curve calculated from the two specific rate 
constants (FIGURE 3). The theoretical curve agreed with the experimental 
points within 1-3 per cent. The average deviation (error) in analysis of 
mixtures of the two esters was 2 per cent (absolute). The experimental 
details will be reported in detail in a subsequent publication. 

Analysis of Mixtures of Olefins. The rates of reaction of perbenzoic acid 
with the following types of olefins, RCH=CH2, RC(CH;)==CHe2, RCH= 
CHR, and RC(CH;)==CHR (mixture of chloroform and benzene as solvent) 
to form epoxides was found to be 1:25:30:400 at 25°. It is thus possible to 
analyze mixtures containing two different types of olefins or other ethylenic 
compounds. As an example of the analysis of mixtures of ethylenic com- 
pounds, mixtures of oleic and undecylenic acid were analyzed.“ These 
compounds contain double bonds of the types, RCH==-CHR and RCH= 
CHe, respectively. The accuracy of the analyses was about 1 per cent 
(absolute). 

The homo-competitive method of analysis was also applied to the analysis 
of synthetic rubbers for “internal” and ‘external’? double bonds. The 
analysis presented several special difficulties but it was found possible to 
overcome these and to devise a practical and convenient procedure. The 
experimental work has been described and discussed in detail in papers 
published elsewhere. 

Analysis of Mixtures of Carbonyl Compounds. Analysis of mixtures of 
carbonyl compounds might be based on rates of formation of aldehyde 
and ketone bisulfite addition products. It was found, however, that it is 
generally preferable to base the analyses on rates of decomposition of the 
addition products. The rates of decomposition of formaldehyde, acetalde- 
hyde, propionaldehyde, benzaldehyde, and acetone bisulfites were found 
to be in the ratio of 1:60:160:1600:400 (at 25° and pH 3.4). As an illus- 
tration of the analytical method, mixtures of acetone and formaldehyde 
were analyzed. A preliminary investigation of the effect of pH on the rates 
of decomposition disclosed that the ratio of rates of decomposition was 
very large in solutions of high acidity. Thus, in 1 N hydrochloric acid at 
25°, the rates of decomposition of acetone bisulfite and formaldehyde bi- 
sulfite are in the ratio of more than 10,000:1. Analyses of mixtures of the 
two carbonyl compounds yielded results accurate to at least 1 per cent. 
This investigation will be described in detail elsewhere. 

It is expected that analysis based on homo-competitive reactions will be 
used chiefly for the following types of mixtures: mixtures of homologous or- 
ganic compounds, especially of methyl and ethyl compounds; mixtures of 
primary and secondary, or of secondary and tertiary organic compounds; 
mixtures of ortho- and para-isomers, efc.; mixtures of different types of 
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ethylenic compounds; mixtures of monochloro- and dichloro-organic com- 
pounds, etc. 


Summary 


A method for the analysis of mixtures of compounds based on the differ- 
ence in rates of reaction of the components of the mixture with a given 
reagent has been described. The method is useful chiefly for the analysis 
of mixtures of two organic compounds that contain the same functional 
group. In general, the procedure requires two titrations for the analysis 
of each mixture. One titration establishes the total molar amount of the 
two compounds, A and B, in the mixture. The second titration gives the 
percentage of A + B reacted with the reagent under standard conditions 
(i.e., with a given initial molar concentration of reagent and A + B, ata 
given temperature, and for a given optimum reaction period). The per- 
centage of A + B reacted (or unreacted) after the optimum reaction period 
is then compared with a calibration curve that gives the original composition 
of the mixture directly. 

Mathematical expressions have been developed that represent the kinetics 
of the reaction of a reagent with a mixture of two compounds under various 
conditions. The computations involved in the theoretical calculation of 
calibration curves are illustrated. (Calibration curves can also be found 
from the experimental results of analyses of mixtures of known composition.) 

Experimental results, of the applications of the method to the analyses of 
several types of mixtures of organic compounds are discussed. 
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AUTOMATIC PAPER CHROMATOGRAPHY 


By Ralph H. Miiller and Doris L. Clegg* 
Department of Chemistry, New York University, New York, N.Y. 


Introduction 


Chromatography is a highly developed subject which provides techniques 
well-nigh indispensable in analysis, preparative organic chemistry, and bio- 
chemistry. The extensive literature is summarized in several mono- 
graphs.” 2-27.28 A special branch of the subject involves the use of paper 
in the form of sheets, discs, strips, or fibres in which very effective separa- 
tions of complex mixtures can be attained. The applications are primarily, 
though not exclusively, analytical, and their value lies mostly in the small 
quantities of material required for examination. The separation of the 
components of a mixture by permitting the mixture to diffuse through a 
paper matrix accomplishes two things. It permits complete separation by 
the simple expedient of cutting out a section of the paper and extracting 
with a suitable solvent, and it permits qualitative identification on the 
basis of the preferential rate of diffusion. Under carefully controlled and 
specified conditions, identification is established in terms of the R; factor, 
which represents the ratio of the distance traversed by the substance to 
that traversed by the solvent. If the process is extended to the use of a 
-second solvent, and in a direction perpendicular to the first, then a new 
and confirmatory set of R¢ values may be used for identification, or to re- 
solve a pair which may have exhibited inappreciable differences in their 
R¢ values in the first separation. This is two-dimensional paper chromatog- 
raphy,> which, among other applications, has been extremely useful in the 
analysis of amino acid mixtures. 

In view of the fact that the techniques of paper chromatography are so 
generally useful, and can be carried out with very simple and inexpensive 
equipment, one may well ask why it is profitable to employ instrumental 
aids in studying paper chromatograms. The use of instrumental aids, 
even fairly complicated ones, are justified for several reasons. For one 
thing, they are capable of refining the present techniques in quantitative 
estimations. What is perhaps more important is the ease with which ex- 
tensive data can be accumulated on the kinetics of the process and on the 
influence of the many factors which affect the process. The simple matter 
of selecting appropriate solvent mixtures is highly empirical and there are 
no guiding generalizations. The busy analyst or biochemist is likely to let 
well enough alone when he has found a suitable solvent mixture for a given 
system, but in all likelihood the next problem will require another series of 
inspired hunches to find the optimum conditions. 

Despite the efforts of several able investigators,®: 4: 6-12: 17-20, 23-26 the theo- 
retical foundations of chromatography are in a very elementary state and 
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practice far outruns the few generalizations which can be used to predict the 
optimum conditions for a new problem. It is likely that ordered progress 
will result from kinetic studies made under carefully controlled conditions 
and that the constant interplay of theory and experiment will eventually 
yield a more logical and coherent description of those complex phenomena. 
This does not overlook the certainty that important and striking applications 
will continue to be made with very simple equipment. 

It is equally certain that the full resources of modern instrumentation can 
be applied to chromatography with profit and that they will be best em- 
ployed if the information is automatically recorded. A series of investiga- 
tions in this laboratory is concerned with the application of various physical 
techniques, optical, electrical, and dielectric, to this problem. 

It is obvious that standard photometric techniques and instruments can 
be employed in the examination of finished paper chromatograms and that 
the results are likely to be more useful and precise as a consequence. A 
recording photometer would seem to offer more complete information and 
eliminate the tedium. Still more useful would be means for following the 
progress of a chromatographic separation. Some experimental approaches 
in this direction are described in this paper. 


Photometry of Paper Chromatograms 


The quantitative evaluation of paper chromatograms has been accom- 
plished with filter photometers, using a special mounting frame and optical 
system for the paper strips.!»? More useful and more readily interpretable 
results are obtained if the measurements involve a transmittancy deter- 
mination rather than a reflectance measurement. This is substantiated 
more fully by the experience of Liebhafsky, reported elsewhere in this 
monograph.* Nevertheless, we have found that rapidly recorded reflec- 
tance measurements have their uses, particularly when it is desired to 
establish the spatial distribution of the material along a strip or disc of paper. 

The difficulty resides in the fact that there is no known simple relationship 
between concentration and reflectivity, whereas, in transmittancy, the devi- 
ations from the Lambert-Beer law are not too serious to prevent a direct 
determination of concentration. In general, it may be said that transmit- 
tancy measurements require greater photometric sensitivity and elegance 
than the reflectance measurements, but are simpler to interpret. We have 
used both and the choice is dictated largely by the type of information which 
is desired. 

A simple arrangement for the evaluation of paper chromatograms by 
reflectance is shown schematically in FIGURE 1. A small synchronous clock- 
motor (1 R.P.M.) is provided with a split mounting arbor on its shaft. 
The specimen holders in this case consisted of selected glass vials upon 
which the paper strips could be mounted by means of Scotch tape. The 
vials could be slipped over the split arbor, where they were held with sufh- 
cient rigidity. The scanning system consists of a 6V concentrated helical 


* Pp. 997-1014. 
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filament lamp operated from storage batteries, a bilateral horizontal slit, S, 
a suitable filter, F, and a stopped-down achromatic lens, L, which focuses a 
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barrier layer photocell, the output of which is fed through a shunt to a 2.5 
mv Brown Electronik recorder. Since the scanning drum and recorder 
chart are synchronously driven, the resulting record is accurately related to 
distance along the chromatogram. The scalar relationship 1 is most readily 
established by scanning a piece of cross-section paper, in which every division 
is registered as it passes the scanning beam of light. Some degree of varia- 
tion in the scale can be accomplished by changing the recorder chart speed, 
but, in this instrument, the scanning rate is fixed. The choice of filter F 
involves a compromise for multi-component colored systems. Obviously, 
it is impossible to choose a spectral region which is best suited to all com- 
ponents present. For a general survey of the specimen, a region is chosen 
for which all components exhibit some degree of absorption. Since the 
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FicureE 3. Recorded paper chromatogram of leaf pigments. 


information is so rapidly recorded and without further attention, this 
limitation is minimized by repeating the scanning operation with a different 
filter. ‘There are some instances in which the filter limitation becomes an 
asset. For example, if a specimen has been examined for the presence of a 
trace of impurity, it is possible to exaggerate or enhance its presence by 
scanning in a spectral region in which its absorption is much greater than 
the bulk constituents. 

Typical examples of recordings made with this equipment are shown in 
FIGURES 2 and 3. In both illustrations, elapsed time reads from right to 
left, the top represents maximum reflectivity and the reflectivity decreases 
toward the bottom. Both illustrations reveal little of interest to the chro- 
matographer, other than the degree of reproducibility, speed of recording, 
and certain obvious defects which were remedied in an improved design. 
In FIGURE 2, each sharp downward spike (56.1 per cent reflectivity) cor- 
responds to the toluidine blue zone. The rise to the immediate left corres- 
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ponds to a clear, more faintly colored zone, followed by the brilliant green 
zone. Beyond this, the rapid rise in reflectivity represents clear paper. 
The irregularities encompassed by the indicator x correspond to the over- 
lapping edges of the paper strip and the Scotch tape binding. This nuisance 
has been eliminated in the improved design by a microswitch arrangement 
which short circuits the recorder for the duration of the overlap interval. 
The elapsed time for a complete scan is one minute and, to show the degree 
of reproducibility, is repeated approximately two and a half times 
in FIGURE 2. 
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FIGURE 4. 


In FIGURE 3, a recording is shown for a chromatogram of an autumn leaf 
extract. The peaks at A correspond to a pink zone of an anthocyanine. 
In duplicate samples, such zones were cut out, extracted, and measured at 
pH 3.6 in a Beckman spectrophotometer, from which the close resemblance 
of the absorption spectrum to reported values for typical anthocyanines 
was demonstrated. There was no interest at the time in further identifica- 
tion or in the stronger absorption shown in FIGURE 3 for three distinct and 
reproducible yellow and brown bands. 

The advantage of scanning in this manner lies in the speed of operation, 
reproducibility, and the fact that the record is repeated over and over again 
without attention. With minor rearrangements in the optical system, the 
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use of a glass or clear plastic drum permits the recording by transmittancy 
rather than reflection. 
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Ficure 5. Recorded paper chromatogram of methyl violet-aniline orange. 


An improved arrangement for reflectance scanning is shown in FIGURE 4. 
The mechanical elements shown have been assembled from the convenient 
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components used by servo-mechanism engineers.* A Bodine synchronous 
motor provides the drive through suitable gear reduction and the sample 
holders consist of a series of stainless steel drums of varying diameters. 
Each of these is fitted with a steel millimeter tape flush with the periphery 
of the drum and against which the strip chromatogram can be aligned. The 
optical scanning unit shown in the right foreground is a commercial device 
widely employed in wire and radio transmission of pictures.f In this unit, 
an area of the sample is flood-lighted by the lamp and lens system. A 
microscope objective with a built-in stop picks up reflected light from an 
area of approximately 0.01 square inch and focuses it on a phototube within 
the case. Two varieties of phototube are available, the conventional kind 
and another adaptable to frequency modulation of the signal. We have 
modified this assembly so that it is possible to use a small barrier layer cell 
as well. In addition, provision has been made to incorporate filters in the 
illumination system. This system and the micro-switch cut-out are not 
shown in FIGURE 4. The improvements in this instrument permit a wide 
choice in sample size and a speed of scanning which can be selected quickly 
by a change in drum size or of gear reduction ratio. A typical photometric 
recording obtained with this instrument is shown in FIGURE 5. 


Dynamic Scanning 


The devices which have been described so far have distinct limitations and 
are used primarily to evaluate finished and developed chromatograms. It 
has been found far more useful to provide means for studying the behavior 
of a system while the separations are in progress. The first of these, which 
was set up in the summer of 1948, established the design factors from which 
a simpler and less expensive instrument can be designed. Since this system 
has been described,'* 4 a very brief mention here will suffice. Advantage 
was taken of the well-known use of confined barriers in restricting the flow 
through the paper ‘to a narrow channel. As shown in FIGURE 6, paraffin 
barriers were embossed on filter paper and the sample was deposited at the 
lower portion of the rectangular channel. An eluting liquid could be added 
in the circular area below it, and, as it ascended, the components could be 
selectively driven upward. At a suitably chosen point, the advancing 
column was scanned by a monochromatic beam of light. The general 
arrangement of scanning and recording equipment is shown in FIGURE 7, 
and FIGURE 8 shows typical results for the preferential rate of movement of 
solvent and components as they passed the scanning spot of light. 

Although our interest and concern has been primarily with dyes and 
pigments, or with specimens from which a colored compound can be de- 
veloped, it is recognized that the recording of colorless substances is of 
much more general interest. These techniques can be adapted to the ultra- 
violet, but in our experience it is not essential, as one might imagine, to use 
expensive quartz optics and phototubes for this region. Some very en- 
couraging results have been obtained with the use of fluorescent targets and 
these may be either in the form of receiving screens or fluorescent patches 


* Available from Servomechanisms, Inc., Mineola, L. I 
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incorporated in a portion of the paper. These investigations are not suffi- 
ciently advanced to report on at this point. 

Means are available for the continuous examination of chromatograms 
during their development. It is conceivable that one can move a beam of 
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light rapidly over the sample and present the resultant transmittance or 
reflectance continuously on the screen of a cathode ray tube. Elaborate 
mechanical means for moving the beam are unnecessary, because one can 
better employ a cathode ray tube as a scanning light source. This technique 
has been described by one of us.!*»!6 The principle of cathode ray scanning 
as applied to paper chromatograms is shown in FIGURE 9. The cathode ray 
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tube at the lower left has its beam deflected in a vertical direction by the 
usual saw-tooth sweep generator. The flying spot of light is focussed on 
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the paper by a lens. The light which is transmitted through the specimen 
is received by a photomultiplier tube and drives the vertical deflection 
plates of a cathode ray oscillograph. An intermediate logarithmic amplifier 
is useful, though not essential, because the presentation will then be in 
terms of extinction values. The horizontal deflector plates of the oscillo- 
graph are driven by the same sweep generator. As a result, the horizontal 
axis represents distance along the paper strip, and the vertical excursions of 
the pattern are a measure of the opacity. The advantages of this technique 
are not immediately apparent. It might be assumed that the cathode ray 
picture should be photographed from time to time if anything other than an 
enlarged view of the process is desired. Actually, the instantaneous values 
of the photo current may be used in several ways. For example, suppose 
it is desired to establish the time at which solvent and two or three com- 
ponents reach one or several reference points on the paper. In this case, 
suitable triggers would be synchronized with the sweep generator and their 
presence would be ignored until the scanned pattern produced fluctuations 
at the same instant (distance). ‘Their simultaneous occurrence would then 
be recorded. No ambiguity would result, because solvent increases the 
transmittancy and colored zones decrease it, and the pulse directions would 
be reversed as a consequence. The electronic means for such temporal 
discrimination are well established and are discussed briefly in reference 16. 

The principal advantage of rapid scanning as described above is the ease 
with which the resultant photocurrent may be modified electronically for 
indication or control purposes. Thus, if the presence of an absorbing zone 
gives rise to a gently rising current, then, at high speeds of scanning, this is 
readily amplified, differentiated, and made to generate an abrupt pulse. If 
the zone becomes more translucent, as it will when pure solvent appears, a 
pulse of the opposite polarity can be generated. Despite these advantages, 
there is some merit and particular use for intermittent scanning. Consider 
the case in which the rise of eluting liquid will occur over a useful range and 
effect measurable separation of several components in a period of three 
hours. If provision is made to scan the strip with a spot of light from 
bottom to top at a uniform rate and within 30 seconds, the high speed re- 
corder can trace the photometric record with adequate fidelity. If this 
process is repeated automatically at 2- or 3-minute intervals, one will have 
a set of 90 or 60 tracings for the entire process. From these data, rather 
complete information can be deduced, such as eluent diffusion rate, relative 
R; values, concentration distribution in the various zones, efc. Since high 
chart speeds are required, it is useful to stop the chart for one or two minute 
intervals and initiate its motion automatically some 5 or 10 seconds before 
each scan. Appropriate markers to indicate the exact beginning of a scan 
are most easily incorporated on the record by means of pencil reference lines 
on the sample. When these are scanned, very sharp index lines appear on 
the chart. 

These arrangements, with minor modifications, are all useful in the study 
of paper chromatography. By their use, interesting regularities and de- 
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grees of predictability can be evolved in what is otherwise a highly empirical, 


but extremely useful, technique of modern analysis. 
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